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When using merely ones and zips 
To carry sound from the lips 
The leading issue always fits 
•IN A SECOND - HOW MANY BITS?' 

That's the question we address 
Lower bit rate is the stress 
The rate depends upon the trade 
Of bits and bucks for service grade 

'M,R, Aaron', in 

Guest editorial, IEEE Trans, Commun. 
vol, C0M-30, No. 4, April 1982. 
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ABSTRACT 

Speech, due to its inherent redundancy, lends itself 
readily to the method of predictive coding. By employing a 
linear predictor and MMSE criterion, the DPCM loop is 
studied with adaptation of quantizer and predictor. The 
autocorrelation coefficients, for updating the linear 
predictor tap coefficients are determined by maximum entropy 
method, a more accurate method for power spectrum evaluation 
than the conventional windowing procedures. 

With variable rate Huffman coding, the quantizer output 
is further compressed in bandwidth. It is found that nearly 
17.9 Kb/s and 12.8 Kb/s are needed as final data rate for 
4 bit and 3 bit quantizers respectively. In addition, 

1-2 Kb/s is needed for parameters transmission from coder to 
decoder. , 

The study is simulated on the DEC- 10 computer of 
I . I • T . f Kanpur . 



CHAPTER 1 


INTRODUCTION 

1*1 HISTORICAL BACKGROUND 

The purpose of speech is effective conununication* 
Nearly 90 percent of the processed signal in communication 
is speech. This enormous figure creates an ever increasii^ 
demand on the existing digital facilities (e.g. communica- 
tion channels and data storage), which can be alleviated to 
great extent by representing the same amount of information 
with fewer bits at the expense of more sophisticated proce- 
ssing. The method by which this is achieved goes by various 
names of ^redundancy reduction’, ’data compression' and 
others. Effectively redundancy could be used to decrease 
bandwidth, to increase speed, to lower probability of error 
or to lower average signal power [l]. 

Historically, redundancy reduction was investigated in 
early 1950’ s by Oliver, Elias and others. Each of these 
papers relied heavily on the theory of prediction developed 
by Wiener in the early 1940' s. According to Wiener [2], any 
operation is to be based on the data available alone, and no 
other additional apriori knowledge is available. That is, 
the prediction is to be based on the past values of input 
only. 




fig* nt BASIC DPCK 
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Speech and picture signal, due to their inherent 
redundancy present, lend themselves readily to the 
reduction techniques of predictive coding* Instead of 
sending the quantized sample directly like in POA, if only 
the variations between successive samples could be sent 
and at the receiver from these variations that are received 
the original highly correlated signail with its redundancy 
could be reconstructed, a reduction in data rate is achieved. 
This is the basic principle of all differential coders. 

Differential pulse code modulation (DPCM) was one of 
the early versions of differential coders to be studied. 
According to NS Jayant ^Simple DPCM is still one of the 
classic triumphs of intuitional waveform coding' [3], A 
DPCM system is shown in Fig. 1,1, 

The predictor, based on past samples, forecasts what 
is the present sample. This predicted part is subtracted 
from the input, and the transmitted signal is the unpredi- 
c table compKjnent, The predictor parameters are determined 
by different methods. Various criteria are employed towards 
making the predictor ideal (i.e.) in making the unpredictable 
component which is prediction error tend towards zero. 

Each method of error minimisation is assigned a 'cost'. 
The standard cost functions are the following i 

i) Q>st function Q(e) ef e is equal to the squared 

C(e) * e^ (Xa) 


error 
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ii) Cost-function C(e) is equal to the modulus of error 


C(e) * |e| 


( 1 . 2 ) 


iii) When over a small range centred around origin C(e) = O 


and in the remaining region equals one 


C{e) 


* 


(1.3) 


The first method of error minimisation goes by the name 
Minimum Mean Squared Error (MMSE) criterion. In this method, 
the estimate of the present sample is always the mean of the 
a posteriori density (conditional mean) [4]. The second 
approach, known as the Maximum a posteriori estimation (MAP) 
gives the estimate value as the one at which the a posteriori 
density is maximum. When the a posteriori density is 
symmetrical about the conditional mean, the optimum estimate 


is the conditional mean. 


The prediction process can be either linear or non- 
linear [&]• Linear prediction is easily amenable to mathe- 
matical formulations and gives closer insight into the , 
problem. Under the linear prediction case, the error at time 
•n' is uncorrelated with the input at everypoint of the 
observation interval. That is they are orthogonal. For 
most of the analysis purposes, the choice for input is 
Gaussian samples. Speech with some approximations can be 
modelled Gaussian. For linear prediction with Gaussian input. 
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MMSE and MAP estimates coincide. For Gaussian assumption 
of input, the optimum linear predictor for minimising 

mean squared error is the best of any type. In other words, 

\ 

a nonlinear processor cannot give an estimate with a smaller 
mean squared error. Because of all these observations, 
in this thesis linear prediction and MMSE criterion are used. 

1.2 PROBLEM FORMULATION 

The DPGM system shown in Fig. 1,1 is simulated on DEC-IO 
computer and studied in this thesis. This involves the 
following t 

i) Generation of a sequence of Gaussian numbers to be used 
as input 

ii) Design of a lowpass digital filter to bandlirait input 
and decoded output 

iii) Determination of the predictor coefficients by auto- 
correlation method (explained in Chapter 2). Maximum 
entropy method is used (explained in Chapter 3) for 
determining the autocorrelation coefficients 

iv) Adaptation of the quantizer to match the changes in 
variance of the input signal 

v) Adaptation of the predictor to match the changes in the 
short term spectrum between adjacent segments due to the 
nonstationarity of speech 

vi) Determination of statistics of probability distribution of 
occupation of different quantizer levels 
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vii) Binary coding of the output of quantizer 

viii) Construction of variable length Huffman codes to 

conserve final data rate. 

The above operations (i) to (viii) pertain to the coder 
side. Corresponding inverse operations are done at the 
receiver to reconstruct the signal back. This necessiates 
the transmission of periodic predictor coefficients and the 
segment standard deviation to the receiver, 

1,3 ORGANIZATION OF THESIS 

After the introductory chapter dealing with historical 
background and problem formulation. Chapter 2 explains 
linear prediction. Also included are the methods to calcu- 
late the solutions and determination of filter stability* 
Chapter 3 deals with maximum entropy method (MEM) and how 
MEM is used to determine the autocorrelation coefficients, 
Basic theoretical explanations pertaining to the statistical 
modelling of speech, optimum quantization, adaptation of 
quantizer and predictor forms the material of Chapter 4, 
Concept of variable length Huffman codes and laie construction 
of the same is covered in the 5th chapter. Conclusions are 
drawn in Chapter 6. Also pointed out are the possible future 
areas of work, pertaining to the thesis problan. Apart from 
the above six chapters are included two appendices A and B, 



7 


Starting with signal flow diagram. Appendix A deals 
with the generation of Gaussian random numbers and the design 
of the digital LPF. This design is explained in a logical 
manner via flowcharts. Appendix B is the DlCFORTRAN-10 
program listing of the simulation experiment. Ample coimnent 
statements are included to facilitate ease of following the 
program. 
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CHAPTER 2 


LINEAR PREDICTION 

2.1 LINEAR PREDICTIVE FILTER 

One of the most powerful speech processing techniques 
is the method of linear prediction. For redundancy removal 
to curb data rate, (bandwidth) predictive coding has emerged 
as a powerful tool. It is applied not only in speech but 
in image processing also. 

The basic idea is to predict what would be the present 
sample based on the knowledge of previpus ’p* san^les where 
1 < p ^ n-1. As mentioned in Section 1.1 we assume a 
linear prediction and MMSE as the criterion. 

The analysis is approached from the principle of 
orthogonality. If x^, X2 ... are the input (speech) 
data samples f an estimate of a random variable y by a 
linear combination is [l] 

y = ® 1^1 ^ 2^2 ... (^. 1 ) 

The mean squared error is. 

p = (2.2) 

where e( * ) is the expectation operation. 

This error is flainimum when the estimation error 




gjg* 2»1 
B^(n) 


Forward jprad 


Predictor of 






If 



8(a) 


ror 
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is orthogonal to the data Xj^i k * ... N [2], i.e. 

e(ex^) = 0 (2.4) 

Then the resulting mean squared error is 

P = e(e 2 ) = E(ey) (2.5) 

Expanding these results to estimate the value of s(n+l) of 
a random signal, we get (assuming the random signal as «s(n)) 

sj^Cn+l) = a^ s(n) + •., + a{j s(n-N+l) (2,6) 

N 

“ \ s(n-k+l) (2,7) 

This involves the N most recent values of s(n-k). 

A 

®n(*') “ - ^(n) (2.8) 

The coefficients a|J that minimize the MS value of the 
prediction error defines an FIR filter as shown in Fig.2,1. 
The superscripts f and N denote the forward error filter 
and its order respectively. e^(n) is called the forward 
prediction error. By equation (2*4) 

r f 

£[ej^j(n) s(n-k)] a 0 1 i k £ N (2.9) 

This yields the system in terms of liie autocorrelation 
coefficients as a set of linear equations 
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r(o)a5 + r(l)a 2 + ... + r(n-l)a{j = r(l) 

r(l)a^ + r(o)a 2 + ••• + r(n-2)a^ = r(2) 

. . ( 2 , 10 ) 

♦ m 

• • 

r(n-l)a^ + r(n- 2 )a 2 ^ 

Replacing (n+1) by n in equation (2.6) 



Taking z transforms, the system transfer function 
e^(z) ^ 

"n(==> ^ - I = c®k *('>->')/ s(n)] 


^ , y oN -k 

as X — £ at. z 

k=:l 


Since s(n-k) = z*"^ s(n). 


(2.14) 

(2.15) 

(2.16) 


If the signal s(n) could be estimated from the recent past 
’N’ values, equally well, it could be estimated from 'N* 
closest future values. This gives rise to the concept of 
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Backward error filter. By a similar approach as in forward 
error filter, we get 

Ji N ». 

b r(%.-k) = r(k) I 1 ^ k N (2.17) 

r^=l \ 

This is identical to the set of equations (2,10). So, 
a*^d the backward predictor of s(n) is 

s^Cn) = s(n+l) + ... + s(n+N) (2.18) 

The superscript ’b' denotes that it refers to the backward 
error filter. The backward error e^(n) is 

N 

e^(n) = s(n) - 2 g(n+k) (2.19) 

Ic— 1 

It is seen that ejq(n) = ej^(n). So, by taking 'z' transforms 
the system transfer function F^(z) is 

N 

H^(z) s= 1 - S a^ 2 *^ (2.20) 

k*l 

It is seen 

hS( 2) = H^(l/z) (2,21) 

The mean squared error denoted by and 

for the forward and backward error filters respectively are 

N 

ss r(o) - S a^ r(k) a (2.22) 
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2.2 AUTOREGRESSIVE MODELLING 

An autoregressive process (AR) of order M is a random 
signal s(n) satisfying the recursion equation 


s(n) - s(n-l) - ... - Cj^ s(n-M) = i)(n) (2.23) 

where ii(n) is stationary white noise with autocorrelation 
and power spectral density as 

6(m) f ^ (2,24) 

It is seen from the definition that s(n) is the output of 
a linear system with input T)(n) and system function T(2) 


T(z) 



(2,25) 


If the system is stable then s(n) is a stationary process 
represented by the convolution 

s(n) =s S h(r) T>(n-r) (2.26) 

r=o 

where h(n) is the causal inverse transform of T(z)', 

This shows that for any k ^ 1, s(n-k) is a linear condDina- 
tion of only the past values of T)(n)* So, 


eCs(n-k) Ti(n)] « O | k 1 1 
Since ti(n) is white noise by assumption. 


(2.27) 
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f 


The predicted value Sj^(n) of s(n) of order N M is 

^n) = s(n-l) + C2 s(n-2) + ••, + Cj^ s(n-M) 

The predictor error e^^Cn) = s(n) - sj^(n) = T)(n) 

This is seen from equation (2.23) also. The resulting 
is e[e^(n)] == . 

The equation (2.28) is of the form of equation (S 


where 


.N 


1 1 k :< M 


0 

,f/ 


M >< k £ N 
N 

Now with ^ ^ 1 V ^N _“k 


ith Hj^(z) defined as 1 - S aj^ z” , it is seen 

k=l 

that T(z) of the AR process is specified as 
T{z) =-4 

Hi5(z) 

So, clearly s(n) is the output of T(z) with input 
If, S(z) is the output power spectrum, then 

s(z) 

H^(z) I^(l/z) 

2 

Since the numerator is o , it can be written 


S(w) == S(e 


- j- 

ll - S e* 
' • k»l T 


(2,28) 

(2.29) 

error 

. 7 ) 

(2*30) 


(2*31) 

(2.32) 

(2*33) 
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It is seen that by predictor, we could remove the predictable 
portion which on subtraction gives out white noise ■n(n). 

The AR process, with input as T)(n) is able to get back s(n)* 
For these, the knowledge of coefficients a^ = c, * k =s 1 2 
is essential which is done in the following section. 

2,3 ESTIMATION OF PREDICTOR COEFFICIENTS 


From equation (2.10) it is seen that the basic necessity 
is to know the tap coefficients ^ , k = 1,,, M = Cj^. MMSE 

criterion is employed to determine them. The total squared 
error is denoted by E [3] = , 

n. 




n=n^ 

n. 


M 


T* r\ 

' “ ■ ‘^k f Cq = -1 I 


n=n 


n. 


= - 2 


n=:n. 


k=o 
M M 

2 .EC. s(n-k) s(n-i) c. 
k=o iso i 


(2.34) 


(2.35) 


(2.36) 


where n^ and n^^ define the index limits over tA^ich error 
minimization occurs* 


Defining 


n. 


f*. 

^ki ^ s(n^k) s(n-i) 


n=n., 


(2.37) 


total squared error JE ±s 
M M 


E = - E 2 


k=rO i:=0 


(2.38) 
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Equation (2,38) shows that the total squared error is in 
quadratic form (i,e,) atmost powers of two in the coefficien 
are obtained. Minimization of E is obtained by setting the 
partial derivative of E with respect to Cj^i k = 1, *,,, M 
to zero and solving. Therefore, 

3e ^ 

= 0 = -2 Cj^ Cj^^ (2.39) 

Since c = -1 
o 

M . 

“ Sk ^ ^ •••» W (2.40) 

The M unknown predictor coefficients are obtained by solving 
this set of M linear simultaneous equations known as normal 
equations. It is seen that data s(n) from (n^-M) to are 
needed. 

Assuming that a sequence of N speech samples [s(n)] == 
[s(n-l), s(n-2) ,,, s(o)] is available, the covariance 
method is defined by setting n^ = M and = M-1 so that the 
error is minimized only over the interval [M,N-1 ]and all 
N samples are used in calculating the covariance matrix 
elements 

The autocorrelation method is defined by setting = -< 
and n^ sa + «> arxi then defining s(n) « 0 for n ^ 0 and 
n i These limits allow as 
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s= S s(n-k) s(n-i) (2,41) 

ns=-«> 

00 

« S s(n) s(n+ |k-ij) (2,42) 

nS3*« 00 

N-l-jk-il 

=! 2 s(n) s(n+lk-il) =r(lk-il) (2,43) 

nsao 

Thus, although error is minimized over an infinite interval, 
equivalent results are obtained by minimizing only over 
[0,N+M-l] as s(n) is truncated to zero for n >. N and * 0 [4], 

The name covariance method is drawn from control theory 
literature where for a zero mean signal defines the kth 
row and ith column element of a covariance matrix. Auto- 
correlation method draws its name from the fact that re- 
duces to the definition of short term autocorrelation at the 
delay Ik-il, It needs to be emphasized that ttiese terms are 
not based upon the standard usage of the covariance function 
as the correlation function with means removed, 

2.4 METHODS TO COMPUTE THE PARAMETERS 

En^loying either covariance or autocorrelation meliiods 
means solving a set of M linear equations in M unknowns^ 
Various standard methods exist for this [5]. Some of them 
are 

Gauss elimination method. 

Grout reduction method, and 
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The autocorrelation matrix, being symmetric and 
Toeplitz in nature lends itself to some special methods of 
solution. The set of equations can be written as 


~r(o) r(l) r(M-l) ~ 


" ^l" 


■r(l)' 

r(l) 


• 

rs 

« 

• • 

r(M-l) r(o) 


j 

o •• 

1 


• 

# 

r(M) 


= A X = B 


(2.44) 


where A is MxM size} X « Mxl sizei and B is Mxl size. 

To solve this matrix equation (2.44) an elegant 
recursive procedure was derived by Levinson [bj. This 
method was then reformulated by Robinson [?]. Levinson’s 
method assumed the column vector ’B* to be a general vector. 
For the autocorrelation case, it is seen that the elements of 
the vector *B* are those elements already present in' A'. 
Making use of this fact, Durbin |8] devised a faster method. 
The recursive procedure of Durbin is given below. 

Error at the start = r{o) {2.45) 

for values of i from 1,2, ... M 

kj_ - Cr(i) + r(i-j)3/e.^j^ (2.46) 



20 


(2.47) 

1 £ j < i-1 (2.48) 

®i “ ®i-l (2.49) 

where the superscript i indicates the iteration number. 

The final coefficients are 

c. = ct^^ 1 1 j £ M (2.50) 

Jf 

The intermediate quantities kj^^ f l£i£ M are known as 
reflection coefficients or partial correlation (PARCOR) 
coefficients as they are the correlation coefficients of the 
branches when the feed forward and feed back filters are 
put in the lattice structure. The minimum total error e^ 
is eyaluated at each step. This is always positive and 
it decreases (or remains the same) with increasing M. 

2.5 FILTER STABILITY 

The feed forward prediction filter is stable if all its 
roots lie inside the unit circle (i.e.) 

•f ■ ^ N —k 

Hk,(z) =s 1 - E a^ z is a Hurwitz polynomial with 
^ k=;l 

all roots z^ lying inside -Uie unit circle. 

When the autocorrelation coefficients are positive definite 
[l] (which is assured when they a:ce calculated from a nonzero 
signal) 'Uie solution of the autocorrelation method gives 
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predictor coefficients that guarantee that all roots 
lie inside the unit circle [3]* But sometimes, round-off 
errors can cause the autocorrelation matrix to be ill- 
conditioned. Then the check for stability becomes essential. 
While recursively evaluating the Durbin’s algorithm, the 
intermediate parameters so calculated give much insight 
into the stability criterion. At every iteration the error 
e^ is checked whether it is positive, i.e. 

^ 0 I 1 < i 1 M (2.51) 

This actually means that the PARCOR coefficients k^ 

1 ^ i :! M must have absolute value less than unity, i.e. 

Ik^l <1 1 5 i i M (2.52) 

This is the necessary and sufficient condition for the 
filter to be stable [4]. 

2.6 OPTIMUM NUMBER OF POLES 

The filter order M determines the number of poles that 
are present in the filter transfer function. The minimum 
value of M that is necessary to arrive at a particular 
required value of prediction error is chosen. This is 
dependent on many factors. In the autocorrelation method, 
it is seen the coefficients of autocorrelation are solved 
in a set of linear equations to get the predictor parameters 
whose number is ± M. If, say after a lag of x, the 
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correlation becomes negligible such that even if one takes 
into account values beyond t and solves, it is found that 
the tap coefficients of the predictor are almost zero. 

For example, for filter order of < M, we may find that 
the filter coefficients almost taper to zero. In this case, 
nothing much is gained by taking an older of filter > Mj^. 

In the context of speech analysis, studies by various 
authors have established that the optimum number of poles 
for the filter is three [9], [lO]. The plot of predictor 
gain versus number of poles, almost saturate after M — 3. 
N.S. Jayant reports the same phenomenon for M 2 and 
he gives results of simulation studies with a first order 
predictor [ll] • 

In this thesis, therefore, the order of the filter is 
maintained at three. The simulation studies are carried out 
for values of M = 1,2, and 3. 
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THE LINEAR FILTER COEFFICIENTS c^,C 2 and C 3 (FOR 3rd ORDER 
FILTER) FOR SPEECH SAMPLES AS INPITT 


Block 

No. 

^1 

^2 

^3 

1 

0.991675 

-0.615314 

0.146366 

2 

0.726429 

-0.237653 

-0.176381 

3 

1.044862 

-0.684106 

0.172282 

4 

1.072947 

-0.666604 

0.061396 

5 

0.918778 

-0.492328 

-0.043336 

6 

0.918002 

-0.389967 

-0.127984 

7 

0.965157 

-0.519015 

0.009055 

8 

1.036632 

-0.416465 

-0.133012 

9 

1.124920 

-0.517800 

-0.046829 

10 

1.185386 

-0.533546 

-0.035707 

11 

1.241805 

-0.556381 

-0.014824 

12 

1.038111 

-0.073998 

-0.276505 

13 

1.249574 

-0.500309 

0.036164 

14 

1.379794 

-0.601899 

0.048331 

15 

1.446032 

-0.664645 

0.059428 

16 

0.328260 

0.164216 

0.114883 

17 

1.120339 

-0.229878 

-0.045775 

18 

0.926917 

-0,044050 

-0.100819 

19 

0.820616 

0.008654 

-0.030265 

20 

0.693041 

0.157544 

-0.102425 
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THE LINEAR FILTER COEFFICIENTS c^fC^ and C 3 (FOR 3 rd ORDER 
FILTER) FOR GAUSSIAN SAMPLES AS INPUT 


Block 

No. 

^1 

C 2 

^3 

1 

0.013515 

0.010682 

0.002916 

2 

-0.089070 

-0.034122 

-0.114144 

3 

0.080943 

-0.115964 

0.258790 

4 

0.124766 

0.192269 

-0.172938 

5 

-0.102693 

0^160536 

-0.154860 

6 

0.135906 

-0.165468 

-0.057347 

7 

-0.121415 

Oi 092721 

0.022726 

8 

-0.306955 

-0.269164 

0.095525 

9 

0.071145 

-0.244118 

-0.110240 

10 

0.213886 

0.132074 

-0.080623 

11 

-0.436016 

-0,173864 

-0.106404 

12 

0.113613 

0,046058 

-0.006265 

13 

• 0.064803 

0.0763.13 

0.047876 

14 

0.054312 

-0.084474 

0,038279 

15 

0.162500 

-0.101040 

-0.002869 

16 

-0;045142 

0.173489 

0.009722 

17 

-0.250265 

0.069462 

0.046227 

18 

-0.037245 

-0.148705 

-0.220061 

19 

0.061558 

-0.181215 

0,220685 

20 

-0.027743 

-0.011237, 

0.141954 
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CHAPTER 3 


MAXIMUM ENTROPY METHOD 

3.1 STATISTICAL DEFINITION OF INFORMATION 

Let a situation in which L different things 1^, i = 1,...L 
might happen be considered. Let the probability of occurrence 
be p.^^ The probability of occurrence of an event is related 
to information. The relationship between the information and 
the probability is 


I = 1^1 in ^ (3.1) 

where I is the information and a constant = 1 when the 
logarithm is to the base 2. 

If the above system is observed for a long time T (T 
being very large), Pj^ T of Ij^ events, P 2 T of I 2 events and so 
on might happen. The total information about the system is 
then 

^total = ^^Pl^ ^ P 2 T + •••> (3.2) 


The average information per unit time is represented as ‘H* 
and is referred to entropy. Thus 


I 

T T to t al 

n = Y" — 


L 

-k, E 
^ i=l 


In 


(3,3) 


Entropy is a measure of the uncertainty described by a set of 
probabilities. The entropy is zero for a system where all the 
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are zero except one, which is unity. That is there is no 
uncertainty. Entropy is a measure of the disorder in a system. 

With this knowledge of entropy, the concept of maximum 
entropy is fitted in for power spectral determination. The 
assumptions made in the conventional methods [l] concerning 
the data that lie outside the interval of interest, constrain 
the data to be either periodic or zero outside (the interval). 
This leads to unfortunate end effects in the spectral determina- 
tion. So, what is needed, is an approach that is first of all 
consistent with the prior knowledge and second, estimates the 
prior probability assignment that describes the prior informa- 
tion without assuming anything beyond. This was elegantly 
pointed out by [2] as the first principle of data reduction, 
which states 'The result of any transformation imposed on the 
experimental data shall incorporate and be consistent with all 
relevant data and be maximally noncommittal with regard to 
unavailable data'. Jaynes [3j then, formulated his principles 
stating that the prior probability assignment that follows the 
first principle of data reduction is, the one that maximises 
entropy. 

3.2 RELATIONSHIP BETWEEN MEM ESTIMATION AND LINEAR PREDICTIVE 

FILTERING 

Most of power spectrum estimation methods depend on the 
knowledge of measured autocorrelation function. With the known 
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autocorrelation values from lags 0 to n, the (n+l)th lag value 
is determined as that which makes the determinant of the auto- 
correlation Toeplitz matrix nonnegative [4]. The entropy which 
relates to randomness or uncertainity is strongly connected 
to this correlation extension procedure. The entropy of a 
band limited Gaussian time series with power spectrum S(f) is 
proportional to 

/ log S(f) df (3.4) 

0 

Under the constraint that S(f) agrees with the (n+1) measured 
value of the autocorrelation function, i.e, 

W 

^ S(f) cos(2-n;f+x6t)df = r{x) | (x =0 to n) (3.5) 

0 

it becomes necessary to find what S(f) has the maximum entropy* 
Burg [5] answers this question by saying that this problem is 
equivalent to that of designing a (n+l) coefficients linear 
prediction filter with the mean square error power ^N+1* 

The linear prediction filter theory is already dealth with 
in Chapter 2, The one to one correspondence between this and 
a AR process is also mentioned. The input to the AR model has 
the power spectrum as determined by equation (2.33), which is 
nothing but the MEM power spectrum. It is seen that -the MEM 
process and the AR process are exactly equivalent. So, the 
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representation of a stochastic process by an AR model is the 
representation that exhibits the maximum entropy. This was 
indicated by [6j, 

3.3 DETERMINATION OF AUTOCORRELATION FUNCTION BY MEM 

Equation (2,23) gives the recursive equation satisfied 
by an AR process of order M, This is 

s(n)-Cj^ s(n-l) - ... - Cj^ s(n-M) = ii(n) (3,6) 

Multiplying this equation (3.6) by s(n-k) and taking expected 
values, we get, for k ^ 0 

r(k) - Cj^ r(k-l) - - Cj^ r(k-M) = 0 (3,7) 

Since e(T)(n) s(n-k)) = 0 for k > 0» 

Substituting k = 1 to M, a set of equations is obt^inec^ . known 
as Yule-Walker equations. The Yule-Walker equations are similar 
to the set of equations (2,10) [7j, 

If k =0 is included into the equation (3,7),. 

r(o) - Cj^ r(-l) - ... - Cjyj r{-M) = e[Ti(n)’s(n)3 (3,8) 

= r(o) - c^ r(l) - - Cj^ r(M) = (3.9) 

This equation (3,9) can be included into the set of equations 
got from equation (3,7) for 1 £ k i M, The final set of 
equations is as follows t 
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r(o) r(l) r(M) 


1 


2 I 

C 1 

r(l) r(o) r(M-l) 


-^1 


0 

• ' • . • 

r(M)v r(M-l) r(o) 


• 

• 

I • 


• 

• 

0 

j ' 


1 

( 

_ 


The autocorrelation matrix is of size (M+l) x (M+1), the 
coefficients column vector of size (M+l), Since, the R.H.S. 
column vector is also of size (M+l), c can be thought of 
as the power output of a M+l point filter. So, equation 

(3.10) becomes 


r(o) 


r(M) 


r(l) 

r(M) 


1 


%+l 

f 

t 

r(M-l) 

♦ 

• 

• 

r(o) 


• • 

♦ 


7 “ 

0 • 

> 


(3.11) 


For a 2 point filter, the power output is P 2 and the set of 
equations become 


"r(o) 

r(l)~ 


~ 1. ~ 


■ P 2 " 

r(l) 

r(o) 


f 

U 

1 

i 


0 


(3.12) 


The prediction filter can be operated equally well running 
backward over a data set (like the forward). This is the 
backward error filter already mentioned in Chapter 2. The power 
output of the filter can be estimated better by averaging the 
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outputs of the individual forward and backward error filters. 
The average power output of a two point forward error filter 
(assuming a data set of N samples), i.e, for M= 1 


f N-1 

^2 - T3-I 2 

n=i 



^ll®n 


2 


(3.13) 


jf 

where P 2 stands for the power output of a 2 point forward 
error filter and is the coefficient to be determined. 
Similarly, the power output of a two point backward error 
filter is 



N-1 

E 

n=l 



*^ll®n+l 


(3.14) 


This is due to the fact that the forward error filter and 
backward error filter are equivalent, having same coefficients, 
The average power of the two point filter denoted by P 2 is 


P2 = i(p| + P2) 

, N-1 

= 2(N-T) ^ ®n+l 


(3.15) 


N-1 


2 n— ju 

^ll®n ' J ®n”^ll®n+l ' ^ 

n=l 

(3.16) 


This power P 2 is minimised with respect to to find out 
^11* ^® done by setting the derivative of P 2 with res- 

pect to Cj^j^ equal to zero. When this is done, the value of 

=11 N-1 


2 S s„ s„. , 
^ , n n+l 
n=l 


'11 ^ N-1 
E 

■ -nssl 




(3*17) 


’n+l‘ 
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The value of r(o) is determined as 

N 


r(o) = w 2 s 
n=l 


From equation (3.12) it is seen 
r(l) = r{o) 

So, the value of r(l) can be determined. Also P 2 is 
as 

P 2 = r(o) - r(l) 

Substituting equation (3.19) into equation (3.20) 

2 

P 2 *= r(o) - ^(°)^ 11^11 ~ r(o)-r(o)c 2 ^j^ 

2 ' 

» r(o) (1-Cj^j^) 

It is seen that 1 

This procedure is extended to a 3 point filter with 
output P^- The set of equations is 


r(o) r(l) r(2) 


1 


^3 

r(l) r(o) r(l) 


(N 

0 

1 

- 

0 

r(2) r(l) r(o) 


**^22 


0 

- J. 


(3.18) 

(3.19) 
calculated 

(3.20) 

(3.21) 

(3.22) . 

power 


(3.23) 
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From the second row of the matrix equation (3.23) 

r(l) ~ C22 3:(o) - C22 r(l) = 0 
From the 2 point error filter it is known that 

r(l) = 

Substituting this into equation (3.24) 

Cii r(o) - C 21 r(o) - <^22^^11 ^ 

or 

^21 ^ ^11 “ ^ 22^11 
The corresponding power output is 

P3 = r(o) - r(l) C2j^ - r(2) C22 

= r(o) - Cj^j^ r(o) C23^ - r(2) C22 
From the third row of equation (3.23), 
r(2) = C2j^ r(l) + C22 ^(o) 

Substituting this value of r(2) and the value of C2j_f 
equation (3.26) into equation (3*28) 


Pod 


' 22 ^ 


The coefficients of the three point filter is then 

[1, Cj^x^/^22 “ 


3.24) 

(3.25) 

(3.26) 

(3.27) 

(3.28) 

(3.29) 
from 

(3.30) 


-^ 22 ^ 


(S.31) 
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Like in the case of two point filter, the forward error power 
of the three point filter is 


N-2 

n=l ' ®n+2“‘^21®n+l ^ 


(3,32) 


The backward error power of the three point error filter is 


b 1 

^3 = l®n-‘=21®n+l - '=22=n+2l 


(3.33) 


The overall average power is then 


(p^ 4 - F^) 


(3.34) 


The only unknown is C 22 * Differentiating equation (3.34) with 
respect to C 22 equating to zero yields the value of C 22 » 

It is seen that the value of C 22 must not exceed unity. 


The procedure can be extended to a four point filter and 
so on. The general equation for the filter coefficient for 
1 £ £ M and 1 £ k £ M-1, is 


and 


^mk ^ ^in-l,k “ '^mm*^m-l,m-k 


m 


= Pm id “ 
m— 1 


c^ ) 

nun' 


(3.35) 


(3.36) 


The unknown autocorrelation values are determined thus by 
recursively calculating the coefficients of the filter. The 
recursive equation for the autocorrelation coefficient is 
given as 
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Fig 3*1 lattice Filtey 
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k 

r{m) = 2 r(m-n) (3.37) 

n=l 

3 . 4 STABILITY 

The combination of a forward error filter and backward 
error filter is shown in Fig, 3^1. This overall system is also 
generally known as Lattice filter. For example, in the figure 
is shown N sections in cascade. Each section has one delay 
element and two multipliers. The coefficients are 

nothing but the reflection or PARCOR coefficients of the 
lattice. It is seen that these coefficients have magnitude 1.1, 
which is a necessary and sufficient condition for the filter 
to be a prediction error filter* By checking the absolute 
value of these reflection coefficients to be < 1 at every 

recursion, stability is assured. This is a consequence of 
Durbin's modification of Levinson's algorithm, as outlined in 
Chapter 2. 

The forward error filter has system function 

hJ(z) = 1 - (3.38) 

for this filter to be stable this polynomial must be Hurwitz 
i.e., all its roots lie inside the unit circle. Corresponding- 
ly, it follows that all roots of the backward error filter with 

H^(z) « I^(l/z) 


(3.39) 
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lie outside the unit circle. A. Papoulis [9] gives a proof by- 
induction to show that the filters are stable, 

3.5 ORDER OF MEM PREDICTION FILTER 


Proper choice of the order of the filter or AR process is 
an important part of the MEM analysis. Various criteria have 
been suggested to determine the same [lO]. It has been 
noticed by [?] that spurious details were introduced into the 
spectrum by excessively large values of order. Small values 
tend to smoothen out the spectrum thus annihilating some 
significant details. However, the final prediction error (FPE) 
criterion of Akaike [ll] is the most widely used. The 
are a set of AR coefficients based on a certain sample realisa- 
tion of a time series. Akaike 's FPE measures the error that is 
incurred on the average when these coefficients are used as the 
model for other realization of the same time series. For zero 
mean processes, the FPE is given by 



(3.40) 


The order of the AR model is chosen as one that minimises the 
FPE, It has been observed to give a good results if order of 
the process is kept limited to N/2. 
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THE NORMALIZED AUTOCORRELATION COEFFICIENTS OF SPEECH SAMPLES 
CALCULATED BY MEM (FOR BLOCKS OF 50 SAMPLES EACH) 


Block 

No. 

r(o) 

r(l) 

r(2) 

r(3) 

1 

1.0 

0.622342 

0.092937 

-0.144060 

2 

1.0 

0.575682 

0.079000 

-0.255806 

3 

1.0 

0.628727 

0.081145 

-0.173048 

4 

1.0 

0.646241 

0.066455 

-0.298088 

5 

1.0 

0.614346 

0.045496 

-0.303995 

6 

1.0 

0.649136 

0.122861 

-0.268339 

7 

1.0 

0.635983 

0.100567- 

-0.223966 

8 

1.0 

0.710651 

0.225693 

-0.195012 

9 

1.0 

0.732754 

0.272176 

-0.120074 

10 

1.0 

0.764918 

0.345863 

-0.033843 

11 

1.0 

0.793901 

0.417720 

0.062192 

12 

1.0 

■ 0.824057 

0.553609 

0.237223 

13 

1.0 

0.847070 

0.588801 

0.348118 

14 

1.0 

0.881156 

0.656502 

0.423801 

15 

1.0 

0.892936 

0.679633 

0.448714 

16 

1.0 

0.442270 

0.360205 

0.305752 

17 

1,0 

0.884131 

0.720177 

0.557825 

18 

1.0 

0.826158 

0.638436 

0.454566 

19 

1.0 

0.808019 

0.647273 

0.507890 

20 

1.0 

0.749619 

0.600230 

0.431588 
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THE NOIMALIZED AUTOCORRELATION COEFFICIENTS OF GAUSSIANS 
SAMPLES CALCULATED BY MEM (FOR BLOCKS OF 50 SAMPLES EACH) 


Block 

No.* „ 

r(o) 

r(l) 

r(2) 

r(3) 

1 

1.0 

0.013693 

0.010907 

0.003210 

2 

1.0 

-0.084264 

-0.016998 

-0.109755 

3 

1.0 

0.049543 

-0.099133 

0.245021 

4 

1.0 

0.114480 

0.186754 

-0.127627 

5 

1.0 

-0.159526 

0.201623 

-0.201175 

6 

1.0 

0.124272 

-0.155705 

-0.099071 

7 

1.0 

-0.131177 

0.105667 

-0.002266 

8 

1.0 

-0.258009 

-0.214613 

,0.230848 

9 

1.0 

0.079090 

-0.247210 

-0.147135 

10 

1.0 

0.231302 

0.162898 

-0.015232 

11 

1.0 

-0.374069 

0.029038 

-0.054028 

12 

1.0 

0.118712 

0.058801 

0.005883 

13 

1.0 

0.074548 

0.084713 

0,059054 

14 

1.0 

0.047254 

-0.080098 

0.029937 

15 

1.0 

0.147789 

-0.077449 

-0,030387 

16 

1.0 

-0.052555 

0.175350 

-0.007311 

17 

1.0 

-0.262832 

0.123090 

-0.002835 

18 

1.0 

-0.004139 

-0.147640 

-0.213946 

19 

1.0 

0.019275 

-0.175775 

0.206372 

20 

1.0 

-0.029485 

-0.014604 

0.142691 
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CHAPTER 4 


SYSTEM DESCRIPTICa^ 

4.1 OPTIMUM QUANTIZATION 

The purpose of a quantizer is to map the input within 
a given range to an output level. This enables a continuous 
range of input to be restricted to a certain number of 
output levels, for onward digital transmission. The 
quantizer is specified by the end points Xj^ of the 'N^^' 
input ranges and-an output level yj^ corresponding to each 
input range. Based on some distortion measure, the Yj^'s for 
a given N^^ are chosen along with the associated Xj^’s. Max[l] 
uses a distortion measure 'D* as the expected value of f(e) 
where f(*) is some function and e is the quantization error. 
If the input amplitude probability density is p(x) then 


D = e[f (Quantizer In - Quantizer out)] 
N. ^i+1 

~ s' / f(x-yj p(x) dx 

i=l X. 


(4.1) 

(4.2) 


where Xj^j_ 
Xi and x.^^. 


oo and Xj^ = - oo and for an input x between 
the output of the quantizer is y . . To minimise 


D, equation (4.2) is partially differentiated with respect to 


Xi*s and yj,'s and setting them equal to zero. This yields as 
p(x) 0. 





(4.3) 
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and 

^j+1 

f'(x-yj) p(x) dx = O, 3=1,2, (4.4) 

3 

The function f(x) to be used must be a good metric, i.e. , 
f(x) must be monotonically nondecreasing, i.e*, 

f ( 0)=0 ( 4 , 5 ) 

and 


f(x) _ f(-x) (4,6) 

If a cost function f(e) proportional to square of error is 
assumed as 


f(e) = e^ 

(4.7) 

Equation (4.3) becomes 


= 0 

(4.8) 

or 


='k = k=2, ..., 

(-4.9) 

Equation (4.4) becomes 


^j+1 


■' (x-y.) p(x) dx = 0, 3 = 1,2, ..., N, 

J 

(4.10) 

That is y^ is the centroid of the area of p(x) 
and 

between Xj 


The above two equations (4,9) and (4.10) are solved using 
numerical methods as they are not readily amenable to any 



45 


standard analytical methods. Assuming p{x) to be Gaussian 
for different Max [l] obtained solutions numerically and 
gave the x^’s and y^‘s in a tabular form. It is seen that 
the criterion used for distortion minimisation is MMSE, 

Speech, although has been approximated as Gaussian 
by various researchers, is not strictly Gaussian. 

Davenport [2] and Richards [3] have made extensive studies 
with a large number of speakers to arrive at a more accurate 
model for the probability density of speech. Following 
their results, McDonald [4j proposed a special form of 
gamma density as 



(4.11) 


2 

where cr is the variance of input x. A slightly poorer 
approximation is the Laplacian density given as 



It is reported that the above two densities serve as very 
good approximations to real speech density. 

Following exactly the procedure of Max, Paez and 
Glisson [5] assumed p(x) to be Laplacian and Gamma and solved 
for Xjj^’s and y^'s. Table 4,1 shows the values for = 2,4,8 
and 16 (i.e.) for binary 1,2,3 and 4 bit quantizers. The 
’minimum distortion' factor of Max was called as 'optimum' 
by them. 
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Table 4.1 

Optimum quantizer for Laplacian density 




2 

4 


8 

16 


i 

Vi 

^i 

"^i ^i 

^i 

Yi 

^i 

^i 

1 

oo 

0.707 

1.102 0.395 

0.504 

0.222 

0.266 

0.126 

2 



CO 1.810 

1.181 

0.785 

0.566 

0.407 

O 




2.285 

1.576 

0.910 

0.726 

4 




OO 

2.994 

1.317 

1.095 

5 






1.821 

1.540 

6 






2.499 

2.103 

7 






3.605 

2.895 

8 






OO 

4.316 


9 


10 
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4.2 GENERAL DIFFERENTIAL QUAI'ITIZER-PREDICTOR 


Speech signal inherently has high correlation between 
adjacent samples. So, the variance of the first difference 
is smaller than the variance of the input itself. That is. 


= 2 <s^> (l-r(l)) (4.14) 

where e (1) is the first difference and n denotes that it 
n 

corresponds to the nth sample} r(l) is the autocorrelation 
between the adjacent samples. If r(l) is ^ 0.5, obviously 
the difference or error has smaller variance. For a given 
fineness of quantization, the quantization error power is 
proportional to the variance of the quantizer input. So, 
by allowing the error to be the input to the quantizer 
rather than the actual input itself, a better SNR could be 
obtained. Or equivalently for same SNR, fewer number of bits 
are enough for quantization of the signal. Generalizing 
equation (4.13) by subtracting a weighted (by Cj^) times 
of the previous sample, 


= <s^> r(l)) 


(4.15) 

(4.16) 


The above equation has a minimum when 
generalization is done by subtracting 


= r(l). 
M 


A further 
®n-lc "n 






49 


and this becomes the quantizer input. This is the basic 
principle of differential coding [6], 

A general block diagram of DPCM coder— decoder is shown 
in Fig. 4.1. Here e^, the differences or error signal is 
the quantizer input, producing e^^ as the output of quantizer. 
The prediction is based on M past quantized values. The 
coefficients of the predictor are Cj^, k = 1, ..., M. From 
the figure 4.1, the quantization error for the nth sample is 



(4.22) 
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To maximize gain the predictor coefficients are 

chosen to minimise the denominator of equation (4.22)* With 
MMSE as criterion, this is done as mentioned in Chapter 2, 

4.3 PERFORMANCE MEASURES 

It is well known in speech analysis that the SNR is 
not a perfect indicator of speech quality and intelligibility* 
This is because perceived quality ■ and intelligibility 
depend on the subjective loudness of the quantization noise 
and not just on the quantization noise power. Due to the 
auditory masking properties of the human ear, the subjective 
loudness is determined by the spectrum of the quantization 
noise' and its relation to the input signal spectrum [7], [8], 
[9]. The SNR does not reflect these considerations. Still, 
it is the most popular performance measure. 

It is common to do the time averaging in eqn. (4.21) 
over the entire period of the utterance or sentence proce- 
ssed, It is not necessary to do soi one can use shorter 
averaging intervals, A recently proposed approach is to 
compute the SNR over many nonoverlapping blocks of data 
within an utterance [lO], The time variation of the resulting 
short term SNR provides an indication of how well the coder 
is performing over various blocks of data* Based on these 
short term SNR values denoted by SNR B^(dB), (jth block SNR 
in dB) a performance measure SNRSEG (segmental SNR) is 


defined as 
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SNRSEG = y r SNR B.(dB) 
^ j=l ^ 


(4.23) 


Infrequent large values of SNR B- tend to show up better 
in SNRSEG than in SNR [ll]. The short term SNR Bj is a 
useful performance measure for coders with adaptation in 
quantizer and predictor in a block fashion. 


The gain of equation (4.22) is also called as SPER in 
predictive systems like DPCM, That is SPER is motivated 
by a desire to evaluate the predictor performance of DPCM and 
due to the closed loop nature of DPCM with the presence of 
quantizer within the loop, SPER is also affected by the 
quantizer. 


Another performance indicator for predictive systems is 
SNRI defined as [l2], [13] 


CEN" 


SNRI = 




M 

S 

k=l 


<s^> 

— ■ ■ ' ■■ ■' - 

■=k " 


I. No. 


BRARY* 

.8.2.1. At. (4.24) 


where the prediction is based on the past inputs rather than 
on the past quantizer outputs like in DPCM. SNRI is inter- 
preted as the amount by which linear prediction can reduce 
the input signal power and hence sometimes used as a measure 
of maximum utility of linear prediction. It is seen SPER and 
SNRI are closely related. Since SPER is computed based on the 
actual quantizer input e^, it is the true ratio of signal power 



52 


to prediction error power in DPCM# However, SNRI upper 
bounds SPER and hence SNRI is useful as annindicator of 
the maximum performance improvement available. 

Interesting relationships between SNR, SPER and SNRI 
have been derived by Gibson [12], 

4.4 ADAPTIVE QUANTIZATION 

The quantizer design has to be compromised between two 
contradictory requirements. To accommodate a high peak to 
peak signal a wide step size is preferred. .At the same 
time to reduce quantization distortion, it is necessary to 
have a smaller step size. Speech depicts a wide fluctuation 
in the instantaneous amplitude level. The amplitude varia- 
tions depend on the speaker, the environment and has a wide 
swing between voiced and unvoiced segments in a given 
utterance. 

To accommodate these amplitude fluctuations, very often 
the properties of the quantizer are made adaptive. There 
are two basic approaches. In the first, means are provided 
to match the step size of the quantizer to the signal 
variance at the quantizer input. In the second method, the 
quantizer remains fixed with no alteration to step size, 
but is preceded by an amplifier with a time varying gain. 
The time varying gain is calculated based on an estimate of 
the short term variance of the block of the signal under 
consideration. Both the methods are equivalent in operation 
[13]. 
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The quantizer adaptation can be either feed forward 
or backward as the case may be. When the adaptation is 
effected from parameters calculated directly from the input 
signal it is known as feed forward adaptation. When the 
adaptation parameters are determined based on the output of 
the quantizer it is backward adaptation. 

Noll [lO] gives an exhaustive account of the adaptation 
with a variable gain amplifier in front of the quantizer. 

The scheme is depicted in Fig, The adaptation parameter 

is to be proportional to the variance of the input to the 
quantizer.. However, since the error signal e(n) is proportional 
to the input, it is reasonable to calculate the adaptation 
parameters from the input signal [l4]. The short term variance 
is defined as 

L 2 

a^(k) = i S s(n) (4.25) 

^ n=l 

where L is the number of samples in the block 'k’. With 

2 

the knowledge of o (k), the gain for the kth block of the 
speech segment is [l4j (G being the gain), 

G(k) =. 5 ^) (4.26) 

That is the gain of the amplifier is inversely proportional to 
the standard deviation of the block under consideration. This 
necessitates the transmission of the adaptation parameter (a) 
to the receiver for every block. This transmission would account 



54 


for a marginal increase in data rate. In the case of backward 
adaptation scheme, no transmission separately of the adaptation 
parameter is necessary. But, it suffers from the disadvantage 
of increased sensitivity to channel errors, since from the 
quantizer output, not only the signal is reconstructed at the 
receiver, but also the adaptation parameter is extracted. So 
channel errors cause an error not only in the quantizer level 
but also in the quantizer step size. 

The method of step size adaptation, either feed forward 
or backward, is also equally popular. Various algorithms have 
been reported [l5], [16], Most of these methods modify the 
step size (for every new input sample) by a factor depending 
on the knowledge of which quantizer slots were occupied by 
the previous samples. In a simple form the scheme operates 
with a one word memory. The output of a uniform quantizer of 
N^bits is written as 

y„=Hn £n (4,27) 

where y^ is the output for the nth sample : 

+ = 1,3,5 2^1 and 

6^ is the step size for that sample. 

Now the step size 6^+1 is chosen to be the previous step size 
multiplied by a time invariant function M of the code word 
magnitude H^, (i.e.) 
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(4.28) 

where the multiplier function M is properly designed, the 
adaptation method serves to match the step size, at every 
sample, to an updated estimate of the signal variance. 

Depending on the update interval of the adaptation, it 
is referred to as either instantaneous or syllabic (approxima- 
tely equal to the duration of a syllable utterance). The 
scheme governed by equations (4,27) and (4.28) is instantaneous. 
Whereas the variable gain adaptation scheme is syllabic with 
the update time varying from 4 ms to nearly 16 ms. 

Employing adaptation of quantizer results in a better 
matching of the quantizer input to the short term signal 
variance. This reflects as an increase in the SNR performance. 
Approximately 4 - 6 dB improvement has been reported in the 
literature. 

4.5 ADAPTIVE PREDICTION 

The spectral properties of speech vary widely from speaker 
to speaker and also within a speech segment especially between 
voiced and unvoiced portions. Due to the nonstationary nature 
of speech, a fixed predictor cannot predict the signal effi- 
ciently at all times [17] . For example, speech is approximately 
periodic during voiced portions. Then a good prediction is 
possible based on the value of the signal one period earlier. 
However, the period of speech varies with time. So, the 
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predictor coefficients must change with the changing period 
of the input signal. This is accomplished by making the 
predictor adaptive. The number of articles available in the 
literature on this topic seems to be limited. Two reasons 
for this are j (i) prediction adaptation is generally felt 
to be too complex and (ii) the performance improvement due to 
adaptive prediction is not substantial [12]. McDonald's widely 
quoted theoretical results says that 60 percent of the 
available improvement from the predictor is obtained with a 
first order predictor [4], But to obtain good quality of 
speech at data rates less than 16 Kb/s predictor performance 
is to be improved by making it adaptive. The most popular 
method of adaptation of the predictor is to update the auto~ 
correlation coefficients and hence the predictor coefficients 
over durations of 4 ms to 16 ms. During this duration local 
stationarity is assumed. 

A complicated predictor adaptation algorithm is presented 
in [l7]. The overall predictor structure is a combination of 
the two stage predictors - one (PRl) to remove periodicity due 
to pitch or long term redundancy, and another (PR2) to remove 
the short term redundancy due to vocal tract characteristics 
and the shape of the glottal wave. The two predictors in 
z-transform notation are represented as follows : 

PRl(z) = Y (4.29) 

M2 

PR2(z) = ^ c. z 

i=l ^ 


t4*30) 
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where y and are determined later. 

The overall PR(z) of the entire predictor is 


PR(z) = [PRl(z) + PR2(z) [l-PRl(z)]] (4.31) 

-M, ^2 . -M, 

= [yz + E c^z“^ (l-yz -^)] (4.32) 

i=l 

The predictor PR2(z) is obtained in the normal way like any 
other linear predictor. The parameters y and are obtained 
as follows : 


M 


1 


is selected to maximize the normalized correlation 



k 

E 

n=l 


[ " 

n=l 


®n ®n-Mj_ 

3 ? 2 , 

®n ^ ®m-M, 3^ 

m=l ^ 


(4.33) 


The optimum value of M^(opt) is found by a search of computed 

and tabulated values of r., , y is set to have the value 
k ^1 


l(opt) 


S ®n ®n-M 

_ n=l 

'opt “T 2 

ni 


(4.34) 


Typically k = 75 and M^(max) is 150. Though performance wise 
this scheme is very good, due to the large amount of computa- 
tions, adaptive predictors of this type are too complex to 
most real time applications. 
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Like in adaptive quantization, feed forward and backward 
adaptive prediction is possible. One of the backward adaptive 
schemes is the gradient method [l8]. It is reported by 
Cummiskey [18] that his scheme, a combination of backward 
adaptive quantization and backward adaptive prediction failed 
miserably. According to him, the quantizer and the predictor 
are significant contributors to this failure. There are also 
later studies which report that the backward adaptive scheme 
is seriously disturbed by transmission errors of the channel[l9] . 

4.6 TRANSMISSION OF PARAMETERS 

The linear predictor coefficients and the quantizer 
adaptation parameter are conveyed to the receiver once in 
every update interval. Earlier methods directly quantized 
these parameters and sent to the receiver. Itakura's scheme 

[20] coded each parameter by 9 bits. It is noted that a 
relatively high accuracy is needed (8-10 bits per coefficient) 

[21] , Typically, for a first order predictor with 9 bits 
coding per coefficient, 18 bits are needed for transmission of 
these parameters during every update duration. If this 
duration is nearly 20 ms, then 900 bits are needed per sec 
for conveying of parameters to the distant receiver. This 
adds up to the data rate necessary for transmission of the 
quantized differences to give the overall data rate. For three 
predictor coefficients 1800 bits are required for transmission 
of parameters. 
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Instead of quantizing the parameter directly and sending 
them, a better choice would be to use the reflection coefficients, 
(though this does not give any reduction in the bandwidth nece- 
ssary for sending these), this is because, positive definite- 
ness of the autocorrelation coefficients is not always ensured 
when the parameters are quantized [22], This leads to the 
instability of the prediction filter* Also, the quality of the 
decoded speech is a function of the 'maximum perceptual error' 
between the decoded and the original speech. A reasonable cri- 
terion is therefore to minimize the maximum perceptual error* 

It is assumed that an accurate representation of the power 
spectrum is necessary for good quality [22], Thus, the criterion 
for optimal quantization of parameters, is to minimize the maxi- 
mum spectral error due to quantization. This is done by trans- 
forming the reflection coefficients to the log area* ratios and 
linearly quantizing them. 

Recently, it is proposed that two kinds of transformations, 
inverse sine and inverse hyperbolic tangent can be used [23], to 
be followed by a uniform quantizer. It is also pointed out that 
the precision with which each reflection coefficient is to be 
encoded varies, the higher order coefficients requiring less 
precision. The bit allocations for each parameters and the 
maximum to minimum variation necessary for encoding the refle- 
ction coefficients are given in [23] . 

*The term originates from the relation of glottal area to the 
reflection coefficients in the modelling to obtain the vocal 
tract area function estimates. 
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BLOCK SNR FIGURES FOR SPEECH SAMPLES 

BLOCK SIZE - 50 •> QUANTIZER - 3 bit, adaptive*, 

PREDICTOR - First order adaptive 


Blocks 

No. 

BlocirSNFTdBl 

1 

13.124343 

2 

13.045203 

3 

15.637971 

4 

16.261671 

5 

15.204613 

6 

14.124889 

7 

16.064677 

8 

16.582332 

9 

16.956900 

10 

15.343028 

11 

14.359510 

12 

16.524025 

13 

16.900963 

14 

16.023976 

15 

16.694060 

16 

9.023916 

17 

16.683720 

18 

16.643762 

19 

16.515311 

20 

16.178700 

Segment SNR is 

15.394679 dB 


61 


BLOCK SNR FIGURES FOR SPEECH SAMPLES 
BLOCK SIZE ~ 50) QUANTIZER - 3 bit adaptive 
PREDICTOR - Second order adaptive 


Block 

No. 

Block SNR (dB) 

1 

15.684679 

2 

13.711535 

3 

15.377128 

4 

16.902858 

5 

16.570192 

6 

15.625226 

7 

15.514239 

8 

16.268007 

9 

16.493048 

10 

15.972446 

11 

16.756858 

12 

17.084521 

13 

17.489089 

14 

17.120869 

15 

16.635226 

16 

8.932692 

17 

16.053062 

18 

17.011274 

19 

16.522627 

20 

Segment SNR is 15.936255 dB 

16.999514 



BLOCK SNR FIGURES FOR SPEECH SAMPLES 
BLOCK SIZE - 50; QUANTIZER - 3 bit, adaptive 
PREDICTOR - Third order, adaptive 


Block 

No. 

Block SNR (dB)'" ^ 

1 

16.644759 

2 

13.969165 

3 

17.372970 

4 

16.878557 

5 

16,144523 

6 

16.382229 

7 

15.597213 

8 

17.487268 

9 

16.475826 

10 

16.229284 

11 

17.262019 

12 

16.640277 

13 

17.115826 

14 

17.131104 

15 

16.899681 

16 

8.812526 

17 

15,784485 

18 

18.177828 

19 

15.915454 

20 

Segment SNR is 16,147992 dB 

16.038848 



63 


BLOCK SNR FIGURE FOR SPEECH SAMPLES 
BLOCK SIZE - 50; QUANTIZER - 4 bit adaptive 
PREDICTOR - First order adaptive 


Block 

No. 

Block SNR (dB) 

1 

22.759317 

2 

21.333798 

3 

19.596347 

4 

20.105386 

5 

18.534467 

6 

18.418346 

7 

18.419426 

8 

21.740647 

9 

20.878594 

10 

20.500616 

11 

20.662686 

12 

21.966126 

13 

21.015103 

14 

22.719909 

15 

22,561036 

16 

14.860628 

17 

22.668761 

18 

22.254577 

19 

21.234538 

20 

Segment SNR is 20.675214 dB 

21.273964 



BLOCK SNR FIGURES FOR SPEECH SMPLES 


BLOCK SIZE - 50? QUANTIZER 

- 4 bit, adaptive 

PREDICTOR -■ Second order, 

adaptive 

Block 

No. 

Block SNR ( dB') 

1 

18.38610 

2 

20.59669 

3 

20.96681 

4 

21.16609 

5 

20.33772 

6 

21.77319 

7 

21.39272 

8 

22,04423 

9 

22.00732 

10 

21.72962 

11 

21.58912 

12 

20.45459 

13 

20.62245 

14 

23.25254 

15 

22.13913 

16 

16.06377 

17 

20.97418 

18 

22.22658 

19 

21.08975 

20 

21,85206 

Segment SNR is 21,03368 dB 




BLOCK SNR FIGURES FOR SPEECH SAMPLES 
BLOCK SIZE - 50 j QUANTIZER - 4 bit adaptive 
PREDICTOR - Third order adaptive 


Block 

No. 

BlocirSNOW} 

1 

18.80518 

2 

22.95931 

3 

20.55761 

4 

20.43231 

5 

19.87662 

6 

19.97059 

7 

21.47560 

8 

21.32536 

9 

21.92312 

10 

23.19771 

11 

21.70337 

12 

21.12495 

13 

21.31313 

14 

21.73176 

15 

22.76065 

16 

15.57234 

17 

21.27013 

18 

21.34159 

19 

21.87370 

20 

22.22883 


Segment SNR is 21.072198 
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BLOCK SNR FIGURES FOR SPEECH SAMPLES 

BLOCK SIZE lOOj QUANTIZER - 4 bit adaptive 

PREDICTOR - Third order adaptive 


Block Block SNR (dB) 

No. 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Segment SNR is 20.777887 dB 


20.787213 

22.028668 

20.841824 

21.927631 

21.439780 

22.239482 

21.769214 

13.081285 

21.969623 

21.694146 
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BLOCK SNR FIGURES FOR GAUSSIAN SAMPLES 
BLOCK SIZE - 50^ QUANTIZER - 3 BIT, ADAPTIVE 
PREDICTOR - First order adaptive 


Block 

No. 

Block SNR (dB) 

1 

14.176498 

2 

15.041225 

3 

14.441862 

4 

13.569500 

5 

12.963721 

6 

14.192600 

7 

14.074719 

8 

14,454479 

9 

14.964948 

10 

14.876275 

11 

15.050324 

12 

13.718763 

13 

13.948370 

14 

13.789433 

15 

12,059241 

16 

12.719600 

17 

14.115815 

18 

14.060949 

19 

15.278959 

20 

12.887230 


Segment SNR is 14.019225 dB 
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BLOCK SNR FIGURES FOR GAUSSIAN SAMPLES 

BLOCK SIZE - 50? QUANTIZER - 3 bit adaptive 

PREDICTOR — Second ordej' adaptive 

tJXOCK 

No. 


Block SNR (dB) 

1 


14.092083 

2 


14.627620 

3 


14.700254 

4 


14.807927 

5 


13.470818 

6 


14.30 5550 

7 


14.073708 

8 


15.526875 

9 


14.251805 

10 


14.971290 

11 


14.448292 

12 


13.839994 

13 


15.113131 

14 


13.457106 

15 


13.200611 

16 


13.737909 

17 


13.874761 

18 


14.120780 

19 


13.696516 

20 


12.939654 

Segment SNR 

is 14.162834 dB 
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BLOCK SNR FIGURES FOR GAUSSIAN SAMPLES 
BLOCK SIZE - 50} QUANTIZER - 3 bit, adaptive 
PREDICTOR - Third order adaptive 


Block 

No, 

Block SNR (dB) 

1 

14.095141 

2 

14.994264 

3 

13.975322 

4 

15.685120 

5 

13.871730 

6 

14.465337 

7 

13.895567 

8 

15.409681 

9 

14.249836 

10 

15.726452 

11 

14.046402 

12 

13.837407 

13 

15.354870 

14 

13.581699 

15 

13.222607 

16 

13.707962 

17 

14.215344 

18 

14.234654 

19 

13.898327 

20 

13.999783 


Segment SNR is 14.323375 dB 
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BLOCK SNR FIGURES FOR GAUSSIAN SAMPLES 
BLOCK SIZE - 50 j QUANTIZER - 4 bit adaptive 
PREDICTOR - First order adaptive 


Block 

No. 

Block SNR (dB) 

1 

19.590659 

2 

19.850695 

3 

20.577953 

4 

19.471216 

5 

18.492211 

6 

18.433614 

7 

19.997118 

8 

19.438855 

9 

18.360900 

10 

19.429696 

11 

20.191251 

12 

18.958665 

13 

19.837039 

14 

19.241767 

15 

19.552083 

16 

19.526083 

17 

18.563282 

18 

20.295842 

19 

18.675315 

20 

19.344851 


Segment SNR is 19.391456 dB 
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BLOCK SNR FIGURES FOR GAUSSIANS SAMPLES 
BLOCK SIZE - 50 j QUANTIZER - 4 bit adaptive 

PREDICTOR - Second order adaptive 


Block 

No. 

Block SNR (dB) 

1 

19.589179 

2 

19.256260 

3 

19.050061 

4 

19.943994 

5 

19.894579 

6 

18.800029 

7 

20.022684 

8 

18.990258 

9 

18,384662 

10 

19.333676 

11 

20.417503 

12 

19.058788 

13 

19.968661 

14 

19.031781 

15 

19.740462 

16 

19.571316 

17 

18,727680 

18 

19.674602 

19 

19.035638 

20 

19.861923 

Segment SNR is 

19.417687 dB 



BLOCK SNR FIGURES FOR GAUSSIAN SAMPLES 
Block size - 50j Quantizer - 4 bit adaptive 
Predictor ~ Third order, adaptive 


Block 

No,. 

~ Block SNR (dB) 

1 

19.466603 

2 

19.140464 

3 

18*220989 

4 

20.079584 

5 

19.384422 

6 

19.630827 

7 

20.697524 

8 

20.411218 

9 

17.702670 

10 

19.637062 

11 

20.307969 

12 

19.498857 

13 

19.443954 

14 

18.249659 

15 

20.015668 

16 

18.785969 

17 

18.990143 

18 

20.148510 

19 

20.783491 

20 

19.592672 

Segment SNR is 

19.509413 dB 
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CHAPTER 5 


VARIABLE LENGTH HUFFMAN CODING 

INTRODUCTION 

The input sample after being processed in the differential 
iback loop, is encoded in the binary format, A normal 
ary code word has a fixed length of n bits/s ample. The 
jth of the code word remains the same irrespective of the 
^ability of occupancy of the quantizer levels. Normally, 
is noticed that the levels of the quantizer around zero 
occupied more often than the levels near the boundaries 
the quantizer. So, it looks logical to use a code such 
t the frequently occuring levels are encoded by a shorter 
gth code word and the less frequent levels by longer lengths, 
s gives an overall average message length which is less than 
t given by a fixed length binary code. That is, for the 
le output SNR, a reduction in final data rate is achieved, 
equivalently, the data rate remaining the same, an increase 
SNR is possible. This of course, is achieved at the expense 
a sophisticated coder with provisions of a buffer to 
msmit at a uniform rate in the channel and facilities to 
ick buffer over flow/under flow conditions. Also to be 
jluded in the coder measures to combat the overflow/underflow 

>blem of the buffer. 
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Historically the variable length code was not actively 

used in communication fields upto early 1950* s, though typical 

example of the variable length code - the Morse Code was 

invented long before. The fixed length code for complete 

decoding needs certain conditions to be fulfilled. If Uj^,U 2 »«.Uj^ 

denote a sequence of L consecutive letters from a discrete 

source selected such that each letter is from an alphabet of 

then the number of different sequences of 

length L is k^. If the code alphabet contains D symbols and 

_N 

of the length of each code word is N, then there are D 
different sequences of code letters available as code words. 

Thus, if a separate code word is to be provided for each source 
sequence (so that complete decoding is possible), then the 
condition to be satisfied is 


N ^ log k 
L - log D 


(5.1) 


So, at least log k/log D code letters per source letter is 
needed for a fixed length code* To quote an example, in 
teletype the source has an alphabet of k = 32 symbols. 
Encoding single source digits (L = 1) into binary code letters 
(0=2) needs N = 5 binary digits per source symbol for 
errorless complete decoding [l]. 

Suppose the k letter alphabet (a^,a 2 , ...» \) has 
probabilities p(ap, p(a 2 ), .... p(a^). Each source letter is 
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to be represented by a code word consisting of a sequence of 
letters from a prescribed code alphabet. Let D denotes the 
different symbols in the code alphabet and denote the 
number of letters in the code word corresponding to aj^. The 
average number of code letters per source letter is n and 

_ k 

n = p(a, ) n. (5,2) 

k^=l ^1 *^1 

From law of large numbers [2], by encoding a very long sequence 
of source letters, the number of code letters per source letter 
will be close to n, with high probability. 

5.2 HUFFMAN CODE 

A constructive procedure for finding an optimum set of 
code words to represent a given set of messages was discovered 
by D.A. Huffman [sj. By 'optimum’, it is meant that no other 
uniquely decodable set of code words has a smaller average 
code word length than the given set. In this thesis, only 
binary codes are considered. 

Let the source letters aj|^,a 2 » •••» aj^ have probabilities 
p(aj_), p(a 2 ), ,.•» p(aj^). It is assumed for simplicity, that 
the letters are ordered so that p(aj^) ^ p(a 2 ) L ••• pCaj^). 
The code words Xj_, X 2 , ..., Xj. corresponding to an optimum code 
and the lengths n^, n 2 , ...» often would not be unique. But 
by imposing certain conditions that are to be satisfied by at- 
least one optimum code, a procedure is arrived upon to construct 
this code. 
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Huffman gave the following as some of the properties of 
optimum variable length codes, 

i) No two code words would consist of an identical arrangement 
of code letters, 

ii) The code words are determined in such a way that no 
additional indication is necessary for decoding, once 
the starting of the sequence of code words is given 

iii) No source letter be coded in such a way that its code word 
is a prefix of any other code word or that many of its 
prefixes are used elsewhere for another source letter* 
Because of this property, this variable lengttx code word 
is also known as prefix condition code, 

iv) If several source letters have equal probabilities (i,e,) 
pCaj^) = p(a 2 ) = p(a 2 ) and so on, it is possible that code 
words determined for these source letters may differ in 
length. But these code words can be interchanged without 
the average length getting altered. 

With the above properties, condition 1 is stated as follows 
For any given source with k L 2 letters, an optimum binary 
code exists in which the two least likely code words Xj^ and 
^k 1 identical lengths and differ in only the last code 

digit. The proof of this is given in [l]. With this condition, 
the problem of constructing an optimum code is reduced to that 
of constructing X^,X 2 » ^k-2 finding the first 

digits of Xj^, Now the source has letters •••» ^Jc-i 
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the probabilities 


p(aj^^) = 1 k-2 

P(^fc-l)+P(a},) k^= fe-1 


(5.3) 


That is the source has reduced ensemble now. The problem of 
finding an optimum code is to find an optimum code for this 
reduced ensemble of source letters. Again, the least two 
probable messages are grouped and a further reduced ensemble 
is generated. Continuing in this way, a point is reached when 
the ensemble has only two messages and then an optimum code is 
clearly designed by assigning ’1' to one message and 'O’ to 
another. 

5.3 CODE CONSTRUCTION AND BUFFER MANAGEMENT 

The quantizer output is the source consisting of 2^ 
different letters. The probability of these letters are to be 
determined. This is done by taking a computer count to deter- 
mine the frequencies of occurrence of different levels for 
3 and 4 bit quantizers, and are given in the results at the 
end of this chapter. With the probabilities known, as explained 
in Section 2, an optimum code is constructed [4], The average 
number of code letters per source letter is calculated. This 
when multiplied by the sampling rate gives the final data rate 
achieved with the optimum code. 
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The variable length code is obviously to be sent at a 
uniform rate on the channel. So, a buffer is required in 
which the samples at the input are of variable rate and 
transmitted on to the channel at uniform rate. This might 
result in situations wherein the buffer might get full or 
totally empty. Both these situations are to be avoided. So, 
an effective buffer management is necessary [5],[6]. 

To combat buffer overflow conditions, an adaptive 
strategy is resorted to. The Huffman coder has two modes of 
operation, (i) buffer normal and (ii) buffer full. In the 
normal mode, the quantizer levels or 2^ letters are encoded 
into the respective Huffman code words. When the buffer is 
full, a buffer status signal causes the quantizer and the 
Huffman coder to switch over to operation with lesser number 
of levels. Typically a quantizer and Huffman coder operating 
with 16 levels under buffer normal conditions, would switch 
to 6 levels under buffer full situation. So, the rest of the 
levels are combined to map into the reduced number of 
levels [5], 

To operate the Huffman coder in the adaptive mode primary 
considerations are the buffer length and the buffer gap. The 
buffer gap is that many number of bits which separates the 
buffer full condition from buffer normal condition. For 
example, let the buffer length be ’r’ bits. When the buffer 
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contents reach this »r» bits, the buffer is full and a 
buffer status signal causes the system to switch to other 
mode of operation with lesser number of levels. This 
reduces the average input to the buffer and causes the 
buffer to gradually empty. When the buffer contents are 
reduced to (r-r^^) bits, the buffer normal status is declared 
and original operation with higher number of levels is 
restored. The quantity r^ bits is what is known as buffer 
gap, an important parameter [7], 

The number of samples in the receive buffer plus those 
in the transmit buffer must always be a constant. They 
represent the delay incurred by employing variable length 
coding. Vifhen the transmit buffer is full, receive buffer is 
empty and vice versa. Provided there are no transmission 
errors, the status of the transmit buffer is always known 
by the receiver. No buffer status signals need be sent. The 
single exception is when the transmit buffer is empty. In 
this case, the buffer must send a dummy signal to mark time 
until it has enough information to send. 
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STATISTICS OF PROBABILITY OF OCCURRENCE OF 3 BIT 
QUANTIZER LEVELS 


No. 

Quantizer 

level 

Prob. of 
occurrence 

Huffman 

code 

1 

-0.222 

0.413 

0 

2 

00.222 

0.299 

10 

3 

-0.785 

0.116 

110 

4 

0.785 

0.114 

1110 

5 

1.576 

0,037 

11110 

6 

2.994 

0.012 

111110 

7 

-1.576 

0.009 

11 11 110 

8 

-2.994 

0.000 

•1111111 


The No. of bits/sample with variable length Huffman code 
is 2.135 
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STASTICS OF PROBABILITY OF OCCURRENCE OF 4 BIT 
QUANTIZER LEVELS 


No* 

Quantizer 

level 

Prob.of 

occurrence 

Huffman 

code 

1 

-0.126 

0.245 

00 

2 

0.126 

0.215 

11 

3 

-0.407 

0.186 

Oil 

4 

0.407 

0.125 

0101 

5 

-0.726 

0.074 

1000 

6 

0.726 

0.067 

1010 

7 

1.095 

0.034 

10010 

8 

1.540 

0.018 

10011 

9 

-1.095 

0.016 

01000 

10 

2.103 

0.007 

101100 

11 

2.895 

0.007 

010010 

12 

-1.540 

0.004 

101110 

13 

4.316 

0.002 

101111 

14 

-2.103 

0.000 

101101 

15 

-2.895 

0.000 

0100111 

16 

-4.316 

0.000 

0100110 

The No 
codes 

. of bits/sample 
is 2.992 

with variable 

rate Huffman 
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CHAPTER 6 


CONCLUSION 

6.1 CONCLUDING REMARKS 

An adaptive predictive coding scheme using maximum 
entropy method is studied in this thesis. The results are 
obtained for a first order, second order and third order 
predictors. Three or four bit quantization is employed. 

The quantizer and the predictor are adaptively exercised. 

It is found that the SNR improvement resulting from quantizer 
adaptation is very nearly 5 dB. This, together with 5-6 dB 
improvement in SNR due to the differential mode operation, 
results in an overall 10-11 dB improvement in SNR over non- 
adaptive PCM, for the same final data rate. The SNR improvement 
due to the adaptive operation of predictor is 1 dB. 

oFor the first order predictor (ADPCM first order), it is 
noticed that the SNR is IQ- 6 dB. This figure is nearly 1,2 dB 
more than the figure quoted by N.S, Jayant [l] for the same 
data rate. The improvement can be attributed to a better 
prediction strategy, resulting from a more accurate determina- 
tion of the autocorrelation coefficient using MEM, For a 
third order predictor the SNR obtained is 21.07 dB for data 
rate of 24 Kb/s which compares favourably with the 21,0 dB 
quoted by P, Noll [2], It is noticed that the saturation value 
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of SNR for order of predictor three and beyond, is the same as 
already quoted in the literature but the rate of growth of SNR 
versus the order of the predictor is fast. This makes the 
present system, fit enough to be studied in hardware, where 
mostly first order predictors are used. 

The concept of ’local stationarity’ is revealed in the 
results. For block from samples 700-800, it is found, for a 
4 bit quantizer with adaptive operation of quantizer and 
predictor, when the block length is 100 samples, that the 
SNRB for the block is 13,08, When the block size is reduced 
to 50 it is found that the SNRB for block of samples 700-750 is 
22,76 dB and for block of samples 750-800 is 15.57 dB, So, 
the block SNR is 19.1 dB, This clearly shows that, as the 
input speech samples show much variations in block of 
700-800, a better prediction by imposing the stationarity 
locally for 50 samples duration is achieved. 

By employing variable length Huffman coding for 4 bit 
quantization, the final data rate is reduced, for the same 
output SNR, from 24 Kb/s to 17.95 Kb/s, resulting in a reduction | 
of the transmission bandwidth needed by 25 percent. Similarly, ; 
for 3 bit quantization the reduction is from 18 Kb/ s to | 

12.81 Kb/s, The bandwidth necessary to transmit the coefficient 
to the receiver, is an additional 1-2 Kb/s. 
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Subjective evaluation has not been done. But it is 
generally accepted that in the evaluation of a system, like 
the one dealt in this thesis, subjective analysis plays a use- 
ful role. 

In the next section, we present some suggestions about 
the possible future avenues of work pertaining to this thesis. 

6.2 SUGGESTIONS FOR FUTURE WORK 

It is already stated that subjective perception of 
speech does not depend on the MS value of the quantization 
error alone. The spectral shaping of quantization noise is 
important. The noise in the formant regions (frequency regions 
where speech energy is concentrated) is masked to a large 
degree by signal itself. So, the frequency components of noise 
around the formant regions can be allowed to have a higher 
energy. This is done by a prefilter in the coder and a post 
filter in the decoder. It is shown by [3] that the spectrum 
of quantizing noise in the reconstructed speech signal is 
inversely proportional to [l-R(w)] ^ where [l-R(w)] is the 
fourier transform of the prefilter impulse response. So, any 
desired spectrum of noise can be achieved to attain a good 
subjective performance. This problem is worth exploring fur- 
ther. The necessary material is found in [ 4 ] , [5] , [6] , [7] . 

The second possible avenue of future study that stems out 
of this thesis is motivated by results of rate distortion theory. 
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The basic idea is that if a slight delay can be tolerated 
at the transmitter, then the speech coder can’ look ahead' 
at several possible paths and hence achieve a better fit to 
the incoming waveform than the single path performed by DPCM, 
This concept of delayed encoding has actually come from the 
results of 'tree' coders [8]. It is shown that [9],[lO] 
tree codes can perform close to the rate distortion bound* 
Recent papers on this topic [ll], [12] promise exciting 
studies in this field. 

The concept of variable frame rate transmission has seen 
rapid developments recently. In the case of feed forward 
adaptation procedures, the parameters are quantized and sent 
to the receiver once in every update interval, also known as 
frame. In the variable frame rate method, parameters are 
transmitted only, when their values have changed sufficiently 
over the interval since their previous transmission. The 
parameters for the un transmitted frames are regenerated at 
the receiver through linear interpolation between the parameters 
of the two adjacent frames. Thus, parameters transmission 
occur more frequently when speech characteristics is changing 
rapidly, as in phoneme transitions and less frequently when 
the speech characteristics are relatively constant. It is 
claimed that variable frame rate transmission can yield sub- 
stantial lower transmission rates [13], The connected litera- 
ture is available in [l4] , [l5] , [l6] . This area deserves 


attention 
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The variable length Huffman coding, though simple and 
elegant has the associated disadvantage of buffer management. 

The buffer overflow/underflow problems are normally difficult 
to handle and involves the use of adaptation as mentioned in 
Chapter 5. A viable alternative is fixed length ’permutation 
code’ [17], It is shown [l8] that optimum amplitude quantization 
and permutation encoding are equivalent in the sense that their 
rate versus distortion performances are identical. Especially 
for small values of distortion, the buffer size in the case of 
variable length code becomes very large resulting in a large 
average coding delay. Though the block length ’n’ of the per- 
mutation code becomes high under this condition, it approaches 
the minimum entropy of the optimum quantizer [19]. The study 
of permutation codes in lieu of variable length Huffman codes 
is also to be looked into. The literature on this subject is 
found in [20], [2l]. 
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Signal flow diagraa 



APPENDIX A 


A.l SIGNAL FLOV/ DESCRIPTION 

The signal flow diagram of the simulation experiment is 
shown in Fig. A.l. The figure comprises the coder and 
decoder sections. The input section of the coder produces 
speech samples or Gaussian samples. These samples are 
processed in the feedback loop comprising the quantizer and 
predictor. Both of them are made adaptive by the adaptation 
strategies explained in Chapter 4, The output of the quantizer 
is then subjected to a counting operation to determine the 
probability distribution statistics of the occupation of 
different quantizer levels. This facilitates the construction 
of the variable length Huffman codes. Binary encoding is done 
on the output of the quantizer and is further encoded in the 
entropy coder. 

Exactly corresponding inverse operations are accomplished 
in the decoder side. The final reconstructed signal is available 
at the output, 

A2,l GAUSSIAN NUMBER GENERATION 

In a computer (like DEC-10 of IIT/K) random number 
sequences are generated easily. The random numbers generated 
have an orderly randomness in them and so they are celled as 
'pseudo random’ sequence of numbers. 
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Gaussian random numbers with mean ji and variance o can 
be generated by initially generating two uniformly distributed 
random numbers rj^ and r 2 » These numbers are then converted 
to the Gaussian distributed random number s by the Box-Muller 
transformation. The Gaussian random number s is 

s = (-2 log^ cos(2'Jtr2) (A.l) 

From the sample s, a sample with specified mean p and variance 
a can be generated as [l] 

X = os + p (A. 2) 

where x is the Gaussian random number with mean p and variance 

o . 

In the simulation experiment, the IMSL routines of the 
DEC-10 system are used to generate r^ and r 2 which are uniformly 
distributed and then equations (A.l) and (A. 2) are applied, 

A2.2 SPEECH SAMPLES 

The speech samples are obtained from 6 KHz sampling. The 
bandwidth is 2940 Hz- The sentence used is 'Have you seen Bill'? 
The digital speech samples are from the Hewlett-Packard's FFT 
analyzer of EE Department, Indian Institute of Science, 

Bangalore. 
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.A. 3 DIGITAL FILTER 

The signal flow diagram shows two kinds of inputs : 
i) speech, ii) Gaussian. The Gaussian samples are to be 
bahdlimited before feeding them to the system input. Also, 
the decoder output is to be filtered to get back the signal. 
These filterings are achieved by a digital low pass filter. 

The filter realization can be done in a number of ways. 

It is called recursive, if the present output sample is a 
function of past outputs as well as past and present inputs. 

A filter is nonrecursive when the present output sample is a 
function of only the past and present inputs. Resorting to 
z transform notation, a digital recursive filter can be written 

N 

2 di • 

-- r - — >‘>0 = 1 

B b.z ^ 
i=o ^ 

where X(z) is the z transform of the input, Y(z) is the z- 
trans former output and H(z) is the z~transforra of system 

function. In the above equation and b^’s are constants 

to be determined. 

The above equation (A. 3) can be rewritten as a difference 

equation by cross multiplying terms as 

N _ ^ . 

Y(z) 2 b-z”^ = X(z) S 

is=o i=o 


{A.4) 
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N ^ N ^ 

2 b.z"^ yY(z) = S a.z“^ X(2) 

i=o i=o ^ 

The term z""^ y(z) is the z“^ transform of y(n-i). 
z ^ of equation (A. 5) leads to 


(A.5) 

Taking 


N N 

^ b. Y(n-i) = ^ a X(n-i) {A.6) 

i=o i=o ^ 

Expanding the LHS of the above equation (A, 6) 

N 

b© + b^ y(n~l) + ,,, = S a. X(n-i) 

i=o ^ 

N N 

Y = ^ a. X(n-i) - ^ b. Y(n-i) 

i=o i=l ^ 

Since b^ = 1. If ^ 1, then 

N N 

'’n = ^ (a./b )X(n-i) - ^ (b./b„)Y{n-i) (A.9) 

i=o i=l 

The structure realized by the above equation (A.9) is called 
the direct form 1 [2] 

A filter with system function H(z) can be realized as a 
cascade of a number of individual lower order filters. The 
H(z) can be written as 

K 

H(z) » a. % H. (z) (A, 10) 

^ i=l ^ 


(A.7) 

(A.8) 
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where each individual Hj^(z) can be a first order section 
like 


H^(2) 




-1 


1 + b^.z 


-1 


or a second order section like 


(A.ll) 


H^Cz) 


1 + a^.z-l t a;.2-2 
1 + b^.z-1 + bj-z- 


(A.12) 


K is known as the order of the overall cascaded filter. The 
realizations of individual H^(z) can be in direct form 1* 

In general the overall H(z) of the cascaded filter can 
be written as 


H(z) 


K 

= 7t 

i=l 



li'^^2i 



+a 

+b 


3i 

3i 



(A*13) 


The filter designed in this thesis is of the above type. This 
requires the knowledge of the following parameters [3]. 


i) Determination of K, the order of the cascaded filter 
(number of individual filters in the cascade) 

ii) The denominator of H(z) 

iii) The numerator of H(z) 

The following are the filter parameters, necessary for 
the design : 

i) - the lower cutoff frequency 2940 Hz 

ii) w ^2 *“ upper cutoff frequency 2990 Hz 
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iii) 

Wg - the 

sampling frequency 

6000 Hz 

iv) 

- the 

ripple in the passband 

0.05 dB 

v) 

p 2 “ the 

ripple in the stopband 

0.01 dB 

vi) 

X - the 

transition ratio 

0.983 


The design of the filter is followed up in flow charts. 
The first flowchart ’FILORD* determines K, This is depicted 
in Fig. A. 2.' Starting with some prescribed values calculated 
from FILORD, the second flowchart *POZE’ determines the 
poles and zeros of the system function. POZE is shown in 
Fig. A. 3. The next flowchart 'FILCON*, calculates the 

coefficients ®2i^ ®3i^ ^li' ^2i ^3i individual 

H^‘s for all i = K, (Fig, A, 4), Also calculated are 

the With the coefficients and K known, the final 

flowchart’ SAMFIL’ (Fig. A, 5) takes the input sample, does 
the filtering and the filtered sample is given at the output. 
To begin with, the initial conditions are assumed to be zero. 
When a sequence is given at the input, it is the flowchart 
SAMFIL which is repeatedly applied as the coefficients are 
calculated once for all at the beginning using the other 
flowcharts. 
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APPENDIX B 


The following pages from 108-138 contain the FORTRAN 
listing of the simulation. The listing starts with the main 
program.% The subroutines then follow. The various routines 
are explained as and when they occur. 



♦ 4 


’r.) 

1 i 


! 1 


Is 
/' a 

r\ r P 
1 1 fl 
t wi tj 

’ ' u e t 
.; !'■ 

, ’1 P t 
. -SPt'. 
■'"i til 
C D 1, 


i-’red ! ct a. ve cn,j.i na«.usl no 'U;;m ai^prifuai. 

IT'"^ aJor.,; A'lth tpelr T/? o-arni'eter s 
* ?Y In the oro.Tri'ii.ln tno 

=-■' thl1o»-ir.Q cholces/seiepfinns ai’^» 

her speech uP oaiussiert sa.j.ciPs are to r»e es 

2 hh a hit cuanr. J S3t ! on is to oe use-^; 
p'p‘ quantiser is to oe hSrte adapti'/s; 
hiiP MPfiictor is tn oe nado ■• 

. hinarv' encooino/cecuiirv'. js 

oinarv codinq, -.vhetner varlrini® 
r-i/LsttcnpUiO i.s to re e'^.oloved; 


h e r 
her 
her 
hpr 


irpn t ; 


aria ot i 78 
to oe owoio^^e.j; 
1 e n 1 1 0 e n t r 0 o y 


'■' 1 >-U\ 

hd r;.r’v 

h,l "-t 

■ Li T iJ A 

f' ( ''!■: 

■..olijr 

'Vl 'r,-'; 

■‘■•AlUrl 

!. 'c 

!d,TO"! 

:)t, 'i,: 

S 1 t.) ! 

iM„ d'. 

''" S I Uw 

: M, ■ • r: 

■A'jT up 

i. hr: 

d T i:r'' 

lU-hi- 

■ oraA 

'■> t'] A L, 

hOfil 

T y:''-,: 

APR 0 


PMU,, 

h \ 

kI ,5. 

■” If' P j.. 

4*:) 

F iO.-i 

'\tC4X 

1 ' *■)/, 

ri\rl\S 

’ r £ 

Yi'JGT 

h''ro(. 

rf5 

F J H >• 

Ai'(4X 

Vi.PiK 

il 

JR 10 

A'!:(4X 

? Fij 

A.-rnsi 

hZCi'. 

Pi-t, ^ 

JA 1 = 

1 0 0 0 

1 -= t u 

') 


J ,uhi.0u) ,Vp(20l 

hpUl you) ,Vidn’?u,2) 


in.j ) ,Oi'. ''I, 1 UO J ,yM'f Ipf iOo ) 

T ^r Vo * I -i Y 5 ^ V } ' ^ t ^ , T 0 a .1^ C 4 1 , T i'?:- f 1 (, y , ) 

'or:.' HT ! I J ;{ t r iT' t?i Is'”' "'■•'■ » '■ “ ' 

o n V ? ! V i ^ v' *: V.' ^ H y 1 H' V M 2 ‘S' PS j ^ o ) , z. r t, r 4 r 2 o j 




» /i I',. M ii j, a r 0 '’i A / u , 0 , i , 0 / 

Si 1 , , 'f R , F.y , ^,c/ij , pc ^ 0 , 0 1 , M , 




f L,.» ' ' n. r.» K* M I r\ .r. u ,i u r\ i j v il H 1/ f 1 

.nTT»''Jh'^ TFrt>Ui‘,!4L,THTS IS DOf’E BY 

rA Ftjh.^n 2iTciTf-j " '*,^5? REOUTHFo.') 


tiT’i' vnAMi’TcjFR fflR ’'4» POp, 4 BTt 


H ^ 4 iU>“4 

0:J i i 1 = 1. p,Q 
''’oJh.'.Tt t j=5,p 


1 „4 

: jRTTu'Jc; 

'f '/ » t;, u 

J 

J ‘4 A {' { i X 


"to,:. 2 7 

> 1 

'■’jRdArC 4X 


’S M 


1 SA-ijiT 

•' ' 

Ca! 
nj . 

'/ 1 '•! ( f , ■ J J * (s , 0 
di*l( i, jjsy.y 

fy f ! 

■p-' 

Rf?; f' u 4 '! t '' r 

c u • ,■ ■■■ i, - i I 


"ifVi PT.LASF KEYih rnK VAf.UF OF iSE'r.E^M 

KFYIfi *’0" png ?nR YEYIM "1" FOR fiAUSRIAW 


' '.• ;:' i { flc , th , a i. , p 2 , ps , k , a , fv 1 1 


v;0 TO 83 

.-/a'CATlftC^rng ' 1 

, i.-'i f' ') 


( 



I ") ') 


'i J * ' ‘i /•' J 'j* ^ » V 1 i'-i , N i , ,.i f. i , ;> T G T ■, , ,. n , I 

" r;;ri -iUt 

Qt'^f .., CL'M ITsI , *r j f f; s ^ «;P5r(, p.,j . 

L 3 '-,(i; sT J =j ,KTL»''='oPt'?Cn,r;}»* ) 


" ' ’ j r-iUt 

<’ i iv ] «> ') 


17 

J 


A T C 

IX 

f 

* 1 ; n 

Yijn 



1 

r 

/ 

162 







»* 

] 

i'? -.1 

A T ( 

4 K 

jf 

*t»' 

YiiU 

hf 


r r 

1 

A 

,1 

C.'A 

'GjT 

r * 

1 / 

r r 

ii 

YJ.'S 

i 

»i , J It 

0 

If 

} 


T 

i 

Pj.; 

i A 

'1 







V 


n.''; 

ATf 


* 

'on 




r 

•r 

{ 


i 6 

0 







} 

,A , 

■' ' 1 , ( 

4X 

f 

' i F 

Yun 


h 

If 

f 

I 

// 

': IT 

iH'.' 

V" I '’ 

tt 0 , «i 

0 

1 



Z 

riv 

'’i* 

,T 

I 7 

FpFy 



/ 


V 

V 

Oi.; 

HO 









.i 

’m ■ 

A r f 

R X 

i 

■' O'"! 

AGJM 

Kf 


1 ^ 

'T 

1 . 

0 ,. 

0 ' 









.,1 

■‘7 

A-i ( 

•iX 

f 

'IK 

YuO 

M 

kF 

1 ^ # 

f 


^.tr 

fi 0 

f! 

f- 

" * T 





.A 

r 

r't 


-’7 * 

, T 

«FC,r4 






f i: 

?5 (•, 1 1 

(:,4 

# 

F ‘v» . P 

) on 

TU 

.jQ 

'riPk: 

0 I 








f;;-* 

.) 

>< '■■. ( 

■V( c 

IX 

f 

'pOG, 

/uO 

Y 

uv 


1 

) 

TC? 

TuT 

'.i G 

r* 

') 






>a',:itgrk td rk /lOuPrTvrv'j 

l'l“, w^iG'Ti?: J!:' TSJr’ijfrj "j 

>. p 1 c r a u i' 1 R f-: rj ^ p j* t >/ p v » 

I H V A ij ' I !•, '} T s ET 1, p ^ ,, 55 » 5 <« 


;TF .-.n'i' 


VjrrcnuT^r;^ VAHIhRlF LFAGI'n EW I'RiiPX , CHU Ti.rG/ 

I' r^'^t \)7 

l'"*f<rTiS‘^”0’' n* '‘"-’li ^AG”K QP' TsiFGFS A5”V' jJf iV'HT 

-.CrtP'.) K.TSFl.t'b 

C J^'T TwijS; 

■iPKi.' ( I'fv: I Ta 1 , FIT. La ' AoC . PUT: ' ) 

2.) 1 lT = ,i,*'lL 

.} •' ■*, * 5 , -J® U I 1 ) ♦ 1 no4 .U , T 1 ,f ^ 1 uO ) 

< » = J-T Li-i )»! GO 

I KTi)=lvGfJ(J) 

S 1 J '. -i V. - ;5 1 ! ;.'j ;,< + !,) 5^ ^ }<; j,^ ) 54 ^ f,; ^ 

J TTPit 

.’ri"in't.sSn;-iM/FLnAT(LJ 
TFf i.SL!,ta ,F.),U)GP 'T,', 

ca(.u AOAPi’(«a«,sPFKic,»'>,.sfc;c;sTu) " 

'"JP'TTftiJL 

IFC I - '^T.DGH to ,1,70 

q A G M ■. : E .„. c V .1 0 W , f'‘ , M A X OR U , AC F C , F PK , G ) 

'.’til iFit K::',,r,),Is1 ,'A,yORD) 

FOPmAK * ftU .gOEFF, ' ,(.'1 2,6/) 

CkhL Sfkm,-, ..•vC..,AC{-f J 

C hLu G'-' ’ ■’ ;■:■ \G' \; , AXCKP , RCF » 

Lt A"‘ ^ c CACF(L,X) ,K = ’ , ■•■■ ) ,!,at , f MAXCftD-l ) ) 

FjR!‘iM'(bF6,4,/) 

CAf.i,. liTNFV '1 -I "r,X) 

|F ( :|SEl.li;2 . '■'■,■■ , ■■ ) „• r.,.T.G7aV 

C p .', .i ..'f- " i. 5-p,:’^.„.Ti'r, n.-:-, 

’ j'; ■ i.’F* r-, I t7.,T3^J,T1 ,l5.r'f.G'5,SLCp^:j, i 

s T;. /' 

t, ,•■('/ ’GO TO 71 ■ ■ ' '■■ . : 




n 

:i ^ 1- f 




71.L 




i 

75 

7 ^ 


J * r r 

. 71 ^; 


r“» 




7 

1 ' 1 7 

77 



J ’ 1 ; 1 

' 1 V 


f"** 

i>h 

I'l' C 

7 '4 


1 ■ 1 ^ 

r. f 4 



7 = 1 c 



1 

7 

1 I 


1 /? 


1 r, 






*•" r(HX, Krij[,J.uWX^n, KP’Si’L.^’S A»fc KDR IH.i ,i' OF ^ tV' iuF. 


' V ■„» •! » Jl .r. f ,l 

'-'rj'fUi ,':>!) 



' , ^ X , 


f J' i 4 ' # ♦ ♦ . 4 : #: ♦ ))t t If ' / ) 
i.riF(i,7j 

F J ■ ' A /A 1 C \i X , ' 'J i'} F U T S A h P i, P ' , 'i X , * 0 ''■» 1' S h I f : F f U’ ' , A 4 , ' 0 . . T 6 P ! • . •' t 
l'Ki-Fl(S> siG;.-H.' ,4X, 'ntCuFFP 
■•Xl'lFttrU) . . 

F ;i P A A ■■’( ( W X , ' * * »■ t ■+• ♦ ♦ ♦ * ' , 4 X , ' ^ * ♦ ♦ » ' , 4 X , ' ♦ >f ♦ jSs t ♦ # X. * » » ' , 4 X ; 

i; = 0,n 

P J ‘v^so , j 

!i:j j.'n;0 I = |,?A 

GU''-' .isSilM’H WGMC I)*»i,CivC r) 
f ?( tSeir-fell .KGjvDGO Tu 1.72 
‘■•'G '*c I JsaGfC IH'StinSTo 

a ,'H Tj, n=PHfu T) 

SA'-rnt J, x)swG>i( I ) 

sa‘-’a'rj. c iJ,l)=uF.cnuT(T) , , , , , 

.." K r PFU U 1. o ) C T , c T ) , yMlI F f I ) , t T ) , KPH( i ) r nfc:CU>!T f 1) ) 

<=")7hA' 1 C13,4x,Fl2.b,4y,F11 , 6 , 4X , F 1 2 . & , 4 X , F I 2 . fi , 4 X , F 3 2 . 0 / ) 

’• ' ■! a 5 i ' ! ( T ) = p G : C 1) - n t C U P t f A ) 
n'i//,\Rsv''‘‘*TTPCl)*’^2 
OV ATOIr^VATuT + OVAh 
thifTK 


1 

C \!.ti iMF'sFJ !,(. A,! 
S4-s!,)nrG.jM,^;n/ 
SiPbiigsio,!)'. ' .... 

-n.tTTFilit^) i/-; 

F jR.-jA'i (r> ;,,**■*■ 

2 7 

:, ; /• 

G 4 T w ® 1 *) * ■'’■'■ >' ’■ ' 
■'fHT.Fl.l jA Oi.n- 

.-A T f 1 ; '' ? 


: 4 L, :/7(; 

1 

j ‘ s' . 4' 

= ’ " ■■4 " “ 1 

4i 



-■.■i ' faiF.-’aiC/f!tJA*gfn) 

s IF' on 

Sm FOP THF AUDVF 
■'■.OATCF'} 

^ jf SUfASlG/pvAX^lf ) 

10 

GAIN OyFfl PC^ IS: * rFl2,6/X) 


^• v'i, r 


CAlaCnT,.,i 


ru • 

•i ' X -•.■ 


FOK .i'Fi;, G 


',<F T;i'"KrF4.u* ) 



i'f ’ffil ,.0* , K'-x, 'Si<R' ) 


v-'M,f>/ f, 


i-‘ I , . « 
r. f 

■■'j'- ■' j f tj f + » + # * + t ', 1 n X ,'*♦♦',/ > 

f.'i' S '■» r,i . . , - A A , ,, 

■ j b‘T T = l,';i, 

• ■'. ' l-'U ,'^i)n ,S 'TU )) 

1 - • Vj r 1 rj ^ , 7 t X r 1 9 ft / ^ 

vJ ’h'''i'~bli.‘i5(aUS.,Ra J 

•’ )'■ ri- He 

n V KL 'U' sS'i - tS Hf-VFl.nATi yj,) 

‘ Hr ) f-H 1 ,fii ) ( 

nvFkr^LF R..P PijR 'fuF F.*'»’iPX I'lrFnVA!, r 
Hi, RSP l"(iU=i ,KTLF='APr.r’nK* ) 

ntM'r.UHlTsI Fil.tE'XaCu.FOH'J 
T rF( 3 , (ii ) ■ 

J VAPinuSDCCUfAMCY HATES HF fRE 0'hV‘fTSF.H*) 

*-KT TF t 1 # 6 / ) 

p 'ul* A 4 h n ^ U ^ ^ ^ ^ ^ J 

F J'^ •tA {-( 30 x ,F 1 2 ,6 , / ) 

CbDSFtllNT 1 = 1 ,t ILFS' APCP.FOP*) 

“/PFi. (llHJTst , FILE-' FHPOR. FOP') 

‘KTl’Fn ,♦)! (0NntTj,,n,T = 1 ,M).1J = 1 ,Mb> 

■.IpFrVtuLlTsl ,FU,E='TwPLT.FUP') 

’ »I = '> r*;) »3 T=t ,ML) 

'TPP:;m(I)M 1T=1 , FILES ^OUTpUL. FOR' ) 

WHnFU ,♦) ((SAMCi{('l(IJ,I),l = i»fO,Tj = i ,ML) 

F IPHAU4X, 'OPEFATL'P TU iMOfErTHE PROGRA^I EXFCRITtJN TS UVEP.FDP 
1 nFTATLFU results UF ThF ArALYlTS PLFPSF HgFEF Tu fHt FTbF 
3 "ABC, FOR AML " AbCu.Fnp" . ' ) 

Sfop 

E,'in ■ ■■■ ' 


Tie niqital flitter routine comorlses of 
\ .Paoroutine DTGFit, 

I I . Function XMAX 

III, siioroutlne digcal 

Tits snorouttne"nif;FrL"is toe filter routine «hicn 
Is used to oanditmit the input, ^ne filter is a fjpF 
vttn tne toilowinq input-outout parameters, 

I jPU'l'S; PI -I s the rtpniP In the pass band. 

R?-is the ripoie In the stoo band, 
rR«is the transition ratio. 

KS-ts the samoHno frequenev. 




1, 1 ^ 


1 b 

?0 


Ot; 


T I 


*’ 'I® flO'i c Hri r t " FTj,ni-^r.«< o^'jlns ftOii 

''.-fi 'h. i IGf-'if CtJ'C , ir i ,^<7 

-U-';'. 5t?b) .n ,'1'V''(2l ) , Z^no ( 20 } , jP ( 20 ) , ( 20 j 


‘ ^ I S f n » 1 0 r ^ r 

= t,0KTt (k1 »=*2) .»/U } 


1 i \ ) 


FF r^sSC-H'^i 1 .-i^2*P2 V f hPoT +Wi ) 
f-‘ K-l=t f.:/.[0 

I 'U3, I 41 3 

r 1 , F P S I , K F i f,. /f'’ U r x , c ' ’ . 

*' J \ , / C!>> 

! = (< 

T = T + I 
';/\-n I 'jssG'i 
'’ypKZ,(0A 

" ■»= 1 1 , F J / t ? . ♦oOh’'' t 1 ) 

Fur2 = ( FS/i .Mi 5'i) F AIA f cniMGA i 

?>'(FPC2.Gr,i- nci ) CiP! tcj 10 
■M = (. 2 + n ) / 2 

’"'le flow chart '’piioiPm’’ ??{v^s nar^. 

1' le f iOA cnart "Fu7.t" Peatns frof^ here, 

Taa folioivlob nortloh heteT'i.lnes ■ the uoi^s aiu ^oros 
ht the transier functjoo. 
y 1 # y £ P S I 

c 1) =SvP{ ( ( 1 .4bA ^i{ 1 > J/2, ) 

O Jhi;.C n=Cf-Pl,X(n, ,X) 

riP!M,7.KP0( 1 ) ,pt>! Lf n 

'■;()''^=CF?i-xcx,o.) 

X = G<M(1) 

F:j?,fc. ( J ) sGiiM^CaOFT ( (<:■•'■ Pi. X,(l . , 0 . ) +P0r.£' ( i ) ) / ( G-iPuX c X , O 
ryf>KX,pui;ti(i) 

TFf hT.Sa.i 3G0 TO 20 
m 15 T=2,?-:i' 

:’.jlGs(,2-»*i)/4 
XsTAMt I) 

TO"-’ = CioPl.XC X,0. ) 

‘Kt 3 5 JsIM'UTG 

7.KROrhG)=S0Pf ( C (1 MXl/Z, ) *( 1 .-Z&FaC J) J / CX-ZEPUC J3 } ) 

ZiiPO f J ) SSURT ( ( ( I MX ) /2 . ) ♦ 1 1 MZFHH ( J J ) / ( Xf iSRh ( J ) ) ) 

PJlt:(r'J)=C^50^.i■U*CMPGX(^^ ,,o, )+Gnx)*(C*'''Pf.'K,f i. ,n, j-Pur pCj) ) 
POr.&Cj)aCSORT( CM ^'Ta!' 

n.'j 25 T = 1 

i}F( jlsRKAf.fpnLF.CJ)) 

Xpr. J )-;"; ' ’.-^'.iiPOLEf J) ) 

Tae ^ m!o X c Oar t " ► "*v:.r.r 

£1.0*' chare"'' ’ M-- 

h- , * . t -4 » .u - 


l + PUtc-f 13 1) 


1 , , O , ; r .J n 


i+3UMJ*(C!4Pl.X( t . ,0. ) tPOUFC JM 


imeis here, , ^ 

" oeglfis from here,ta.ls incXiMes the 


.'h' 

t,> i ! •-' ' 

» -f - 1 
- j i. , 


portloh determines the filter cooffs, 


I 



u. 


■.) J t -?.hPlJ f i ) i: t2 ) 

■ |^=/V(i,l ,1 3 

0 jc tpt, r j *4^7 txp( 

a( T , ;,1 j=iin+K:r;4.i-.{li. 

=}n , 2,2 3=2. *(un-r,uM) 

t f , 2 , ,3 ) =i.) fi >" ?■!{•; + Qfi i4 

n. 3 7 « A ( I , ? , 1 3 + A ( 1 , 2 , 2 3 + A c T , 2 , ? 3 3 

P'jR'w^^l ( Tj ',Kj2.6) 

C )-'■('!. pit, 

r = T 

Xi<:=x '■’/,)( CK, (Jx) 

1=:,Tf1 

TF(f„( J3 .Ky.XF3 GO TG 40 
0:i XD X5 

f-H 1 )=K(J )>!=( ] .-»i ) 

ofc;*,h:cJ) 

F .j 1 1 


F’JMCTIO'J X-'^aX(C,T.J 
lUMKMSr-JM 7C2U) 

K''!.\X=0,iJ 
0 3 50 ! = .1.l 

f ruif3S(xMAX3-AHs(:c(x3).GF.n,ojGL) m 5u 
KiAXsAHSCCCi)) 

TXPt; 49.X«»iAX 
FJRMATf ‘XgAX* ,F].?,A J 
CJ''!.l’ruU[;. 

‘’riTURw 
F '3 f* 

'1’rve €iow chart"FTLCON'*ef3ds nere. 


•I'ae tLOH chart ''SamfTi.'* starts trow here,. 
‘jilBROOTl DTgCaI. I A , E , ViK , rrf , MAd , SIGHn WG . 

At2u,2,3) ,E(2y) ,VIM(20,2),Sir,If^(1 aou) 

tj:pe 5 3,v 

E;)RhAT(*V is', FI?. 6) 
r>iD 65 .1 = 1 , mah 
^ fsisicaMCj) 

OU 60 1*1, NT 
rtPF*,AU.2,l 3 
PJ sV/A(i,2,i) , 

r ■{ f' «.' •»• , p i , p V p 'i 

f’= C ( //. C i , 2 , n 3 •( A( T , 2 3 / A f :i. , 2 , 3 3 3 *V Til ( t , I ) 

! - (aa 1 , 2 , X 3 / A t f , 2 , n 3 f ’ ’ I ? ( T , 2 ) 


,wGfrftion()) 



; ac 1 C 1 -:' 1 : K 1 i i- 1 i: i? 1 . 1 


FjWlM'C’^o Ts',Fl 7 . 6 ) 

u M r ^ U f, 1 ' H * ^ t ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 

r! ;:^i t,2) ) 

t'j': rT.'jne 

TCO|;:*,SiGiM(,7j,.ir,iv(,n 

7 J '-f 1 I'i'lc,. 

Hn; r*jRs 


4 0 


AHAPTf ^--Grvi ,,S!r(V|S ir, , , SEGSTDI 

1 , *i' u 1 d; t fitnii sd'T’ples oefure a^antlan 

SuMSTG-u^rjoies LHe cotai enerqv In t*^e soat^ent 
•■i-Mo,u€ 1 nput , sanc) J er; ta the sp-ansnc 

., StiGiiTi;'- Trie stanHard deviation of tae sednent 

• «i’FfT'>s '-'iGrt-UntPut sa'^oles after triru* tne anolitler witti 
,, , , gain of I/UkGSTo) 

'U '"’Kisr dill 

tT, f, r u S s 0 H T (. s u ' ■ ns T G ) 

nil /p r = 1 , M 

iiartCdCl )/SKGSTi) 

’’ijSTThnti; 

PiiTURiV 

GdO 


1 i5 


s i j s K n u T I d K r F ‘ n-:' ^* ( « g , m a x □ r o , a c f c , pp e , a j 
1 -.'pn rS! '•'Gs-Tne innut samples to cai.the autocorrelation coeffs 
H-The ao.ot samoies in tne segment 

"lAXOnn-Tne max.opder of the Prediction filter of tne 
MFP aloorlthm 

h'JTpiJTSi ACFC-The calculated autocorrelation coetfs 
FPE-AKa j !ce ' s tlriai prediction error 
G-The filter coeffs. of the nredictlon filter 
nihfcpJsrOU WG'U3uOy),ACFCf250),FPF(2S0),GC250),P£R(?5h) 
hlMFMSTUh PEF(250) ,Hf250) 

S;.i ', = 0 0 

DO lib isl,M 

S a |‘1*SILJ M + WG M (. I ) ♦ wGm ( I ) 

coTirivOE 

AeFC(l)«SUH/FLOATlM) 

p-^aAcrcaJ 

OMsacrr ( I ) 

^PP(. i i = FLU ” i- (H + 1) /FLOAT (M-U*PM 
■••F” C I ) = ii.u '■ 



! i;< 

1 'h! (' 


! 


t jO 

I -. 1*1 
1 4 ^’ 

1 'j ’ ■* 

1 


1 v,)''» 

« 

m 

? n V) 

r> 

C 1 ^ :* 
1 ul 


n-'f ^..■■'"■.o) gh t,j i?u 
^).) I la .1 = 1 , iv; 
f’r:'' l,I.)=0.U 

*' J rlTwUe 
n-i 4 ,i=:G.U 
■S ()€;-’ =0 . U 

T J = -, 

'VJ 4 25 J = 1.J J 

V. n 4S6 )U '1-2 . 0 ( v.Ga f vj + .v + 1 1 +PeR ( ,T 1 ) ♦ f wGG f j U Pr.r? 
G;*-: ..sbOGr.H ( f j + „;f i ) tPerc T j 1 + t '‘g?u J j tp'CKf 
.. j V j' riGit-i 
'1t‘'!:'4)=s:jn,4/,sner.^ 
rf f cn TO 140 

0 3 IJ'O J = 1,^J 

'UJft i=u( j^-n+GtN'iUCtK j 

c jH'tont; 

03 135 J = l,,v 
C,i H\ ) = H(J+l ) 

Cji’MTrtllc; 

TJ=JJ-1 

orj 15U J = l,.M 

R£R ( .1 ) = P!=:k C J I +G c % ) ♦Ph.R ( J) 4G C t , '! ) * ..G'l f J ♦ ni 3 - 1 j 
PfiFt JJ = PEF( J + 1 )+(URwl*PfcP (J + 1 l + Gfr-'O^^GOl J + 1 ) 
r:j'’i'T'.!iiK 
,$ ll M s 0 . C! 

OJ 155 J =?,W^J 
S J '• 1 = S U : 4 - A C F C f 3 W + 1 - J ) * G ( J ) 

G'JFXlf'itlt 

ACFCCoHlsSIP'^ 

0^1 = (l ,0-GC*G4)^*;^)+Pi,i 

P .4 = ilM 

TFf^^.^l.Ei^,0) GO TO 100 

FPFCf''H)=FLOAT(M+NtO/FLnAT(V-\'i 4 )*Py, 

FpEtMNl»FPE(MN)/FTE0P 
FPF C Hfv, ) = AbOGI U (1 Pb f NW) ) 

TlfPE*, ACFC(WM) 

r O'^'T iriHi: 

n-J 160 J='1AX0RD, (MAXuPU+l) 

SU"^ = 0,0 

r>j 170 1 = 2, (MAX.OPO-l) 

5iPVi = Sl03-ACFC( J + l-'II^GlI) 

C J’^'TTxPE 
ACFCCJlsSlF^ 

CO^TTfiUE 

no 101 l=2,tMAXORtf 1^ 
ACFCCj,?sACFC(l)/ACFC(l) 

COMTTMIE 
ACFCfDsI .0 


E 4 D 


1 1 5 


( 1)3 

j) 2 



1 ’ D 


Aiiv I ACFC ,.'-AAnKn , acf j 

- .'■ >1 1 K.-fhe au^ocorreel atioa coefts frnTi the i'AF,4 alaorlth’? 
n ,r, f'ax.order ot the autocorrelation luntrix 

n autocorrelation catrlx In tne Toeolitz fort” 

nin-ZlMu'' ACFCC250) ,ACFC7S0,25')) 

' ue routine ojts a ro* vector ( ACFC f 1) ) into 
an nrravC AntCL,,-!) j. 
r, = 1 

fC ^ 

i= 

25o 1 = 1,., 

:-l .', = K-] 

0 = ACFC(3 ) 

It(r^.i-:o.2)ni;i jn ?io 
I- I Tfi 2 15 

2 f i'i 

A^FCtJj 


:2, {K-1 ) 

:ACFC(1T J 


7 If 

nj 22f 
'' t’P(f.,.U = 


Ca-XJi f>£ 

r,:) TO 7 HO 

2 1 h 

r.' /j E F *• 2 

DJ 220 IT 
ACP'lF/'iO 

M U ™ S'J 0 •• 1 

22h 

FJ TT.VitF 

'^■i ] = K 

03 230 il 
A : F a. , 2 3 ) 
= + 1 

2 3D 

njuTTfint 

2on 

li-h+ 1 

KsK + l 

25D 

CjT'TTi'iUF' 

JJ = 2 

f)J 290 1 = 
(M+ 

ij=j jf 1 

2yu 

CJUTINUF 

Rfc:TUR ,4 


:ACFCC1T) 


1 5')=ACFC(JJ) 


FID 


.5UBBnuTIM£ LlNFO^ ( ACF , MAVuRD , K) 

Tfils routine tlnds the solutions of a system of linear 
e..|iiatlons by the (lauss elimination metood . - 

t'lPUTs: ACF-fhe autocorrelation matrix in the roeolJtr. form 
MAXflHD"Tne max.order of the autocorrelation matrix 
nuTpUTSiX-Tne vector X contains the solution of tne system 
linear equations* (i .e) the values of the ta,> -fain 
of the feedoacKf linear oreiictlon filter 

Q |l_ I j ^ ft ^ iW'i, \ ? j(__, _ . ' 

n i4''’5bTaA AC*^(250,250) ,eCF'(25v,25a) ,X(?50) 

') ,) i'h.’- A= 1 # 1 ’•‘“I ) 


O k/1 



1 1 1 


•5 / l i 


IJS 


^.111 

^ i? 
1 ,:j -1 

M'S 

m 

4 

^b<'i 


■^t>7 

15*; 

1(>>! 

1 7 '■' 


^lX=\DStACF(K,K3) 

' ,j ’■1= I K 4 1 3 , !•, 

‘ ’ I I ' ’ ^ 5 

;;^iX = AbS(ACFt».l,f<) ) 

■ ri;;iic, 

Trf Z.r,O.K)Gi; IC* 13^ 

'"J 51v' L=K,(, G + n 
!'i4‘'b = ACFCK,i,) ■ 

" |•■TT.^'!ii 
C j''’ TfrillK 

5.) i3n istKtl),.< 

7=ACFl t,Fj/ACF(K,h3 
')J 3 111 J = k,,CN + 1) 

C:,ti'iTTi:HF 

= .0 

1^2 A=l,iM 
'*ljr,,=:!*'uf.*ACF( K ,KJ 
rj-'ifi ftut: 

V’ I '^‘-' -1^’*' 

K T 'i f:, 1 1 *344) 

F jP'iA r( M'iiF TRlANG'fLAH ‘■iAT«'UX T5*,/) 

f , J=) r'-lAAORD) *1 = 1, C -UKURD-l 1 ) 

FjRMATCbFB.b//) 

True ftacit suBstlttJtlon starts i\ere, 

XlMjeACFCrt, (rHl J )/ACF(M,f' j 
7RITF(1 .34§)r.i,Xt!ii) 

FjRy.ATC ‘the root X' ,Ti, 'TS*,F12.6,/) 

>J i5b 1=1, (M-1) 

»rsM-T 
s tj *-1 SB n , 0 

RJ 350 J=CK+t),N 
•SUMsSUfHACFtK, J)»XC 0 ) 

C j."l' I ARE 

XCK JstACFCK, Cw + 1 ))-SUiv!)/ACF'(K,K J 
''hTTECI ,35?)K,X(E) 

FJR^ATC‘THF root X' ,T1 , 'TS* *El2.b,/) 
n)M:rir4!je 
OJ ro 370 
•vKTTEd .365) 

FJ«MAT(*THP given matrix is SINGOLAH') 

CONTINUE 
R'CTuOA 
K 4 0 


'3.)iO.:niJTl'?E FEOBCKCM, MAX0 RD,;s^GM,X, ERR, QM i’,OHtTR, COUNT, IO uT 
1 , 1 1 RANS, lEwtRO , IFTTRA , MNl , 0 'TSlC. , I SELE4 , 1 SSLER , TEGSTi), 1SELE1 i 
Tills routine can op Rescrihed as the heart of the system 
rr. linearly prectlcts the present value of tne sample oaseo on 



I 1 d 


SUl 

Su ^ 
s 1,1 




s 1 n 


‘it'? 


52 *? 


oast outnuts. 

r ii'if's* H« I rio , of ttiP sdf^ples in L‘>e se j'^ient 

MAKilan- ipie it aenntPs ttie no oi oast output sarriOieso 
t'^AXuRU-Dover ^fnlch tap oredtctlon is based 
•■a'iii-The input satiples of tne seoneot - 
X-'foe linear rreiJctioa filter tao oains 
v^MTSiZ- Size ot the ;^uantiser;3 or 1 hit 
-Se ,1 ect 1 on for nuaatiser ajaotatioa 
istitKi-Selectlon for oinarv encouln i/decoainq 
l.S(j;r.ii-:5-Se1 ertlor for variaole ienoto cohinvi/jecotilafi 
SEuSfh-Tne standard deviation of toe seoiienr 
't.ip aoove 3 paraneters are iceyea in by coe ooeraror, 
o.jTp I'i s : FKR-Tne ouantisatinn error 

O.tT-fne output of toe quantiser 
OMTiP-Tne IfiPut to the unantlser 

CufJUT-Tne vector contalnlna tae quantiser occunancv 
statistics for tne entire Interval 
TiJUt-is the binary encoded vector ot the quantiser 
output 

Ii:PAMS-l.niiT expressed in 2 dimension for the segment 
under cons iderat loo ^ ^ 

lEdiRu-ls the variaole lenoth coded vector of tne 
binary code 

IFTTRA-iFhTHn expressed la l dlnenslon tor the segment 
under consideration 

_ ^ ^ ^ - - .1 *1 s c% /•, nr ^ jp| ' n* m jEk 1/ O r* t" n T* 

VfUTTpr 10 
PAf iOJ, 


under consioerarion 

dgi-Size of tne variaole lenntn vector 

XtnBKllOO) ,£RR(10U) ,WG^n30hu) ,X(?Q),0dTllU0 3,V*^ 

■U •it', ISI U'i COUivTllb} , TOhif 4) , 1 rRA'''S( 100 , 1 ) , Te.''' f Pu( 1 ) ,Tl' 11® 

•’I i 1.^ hi c f n I . n r f 1 1 i 1 


:)ii''fc'‘siSTQM ALOCClOu) 
r i rbcieH Ohisiz. g^iXFui 

■ !■:» A fC» ' . f'k \ r'* n m n 


If ClSELfc:4,Eg,o5Gn TO 50/ 

FOhJiA I’Uox.! UsB' ,4X, 'HSB-l *,4X, ' 1 ' # 4 X , ' bSB ' j 3 X , ^ BIT- i % 4X , 

1 *!n't-2 % 4X. *Bll-i' ,4K, 'BiT-4' ,4X, , 4X, 'BlT-6 * , 4X, BT I - / 

r!,rn'TN’ifc , „ 

TfdShbtl .EQ.UlSEGSTDcl .0 

hj 500 l«l,f4 
IFd.EO.UGO TO 505 
GO TO 610 

0-JT1Ms'4GNC 1 )/SEGSID 

o.xriPcnsafiTiB 
GJ TO 525 
CiprrtA-ue 

FJBKSU.O , ,, 

nj 515 J»1 » t VA t*;Ri}'l ) 

Tf(tI-J>jEO,0)GQ TO 520 

XF))aKlJ)*Ab5C(i:-J) 

FaRKsrOBK+XFO 

CjOTIhUE 

Ca*>ifXNlJE 

Xr0BK(T)=/0!:?F 

f f --sir 

B\CQdi6iZ,EV.3)Gg Tu 550 


’ V.h i'n\SiV4 ( OhTl^^ , OilTUrnO 
'Abu eouM-HioMTaWT^cnufiT) 



1 1 ^ 


f isk!,!*:4 tl- 

: \M i:.’’Cni)4 ( 0 -iTvlMT, iPi!'’’, 1 t'jHI-O.) 

I T,i b?o ■ 

. , . . 1 , r- r- ^ <8 r* , . i 5 ^ - ; *»* *• w 


i > 


B / 1 

B ,< 

B 

B ) 


BS'' 


B ih 

B4''> 
B b B 


S S 4 


S47 

Sue 


B 3 ^ 
B <7 


B it 

Btv',) 


^^Kjrl'tl ,5'^7)?XOuT(I.) ,i,ri,.;*..?isizj , ( T tiM i R J ( J ) , T=1 

<; IP ;yi f 1 u, I I ,7x,I1 ,8 a,1} ,7a,X 1 ,nx , n ,riX,Tl ,«T, Tl ,8X,Ti,dY,T 1 
I c , T i , H X , T I , / ) 
i.n BSB 

r i’i'i' I uMr; 

! ,B30) ( inuTxi,) L=i , grt fSl?.} 
f t ,7X,J1 1 ,7 a, ii ./} 

? '\ r . i., V> " i 5 a 3 ( 0 ; : T X "1 , 0 N T Ll n T X 
7A,f,t, C0Uf’'X1(QNTt)TPr,cnu''>r) 
t F ( i -St?Jih4 , Ey , U ) Tij bBb 
Ca'G, EMCODBCO^IT’u'tT, i<BuT, IEHX’Rj, V'! t 1 
r F r I SGGtiS , EG . 0 ) GO -f-G b 3o 

JBX i'Etl ,557) ( lOi.iTd,) .G=l ,GMi'SiZ) , XTEMXRDX J) , T=t , M O J 
r IPOTClRArll ,7x,i1 r«X,ll ,1 lX,n,BX,Il,dX,Tx,BX,Ti,«X,Tl,3X, 

1 It ,BA,X1 ,/) 
i‘‘,t 5 55 
".iTi r^'ic; 

■J . { r T F. ( t , B i b ) c 1 ri u T ( T , ) , = 1 , 0 '71' S 1 Z J 
rMP iM'(l9A,Xl ,7X,i1 ,»X,Xl ,/) 

7 J-' I FnitK 

<7 JTj. l)=:G’''JtDUT 
F:<tu r)-Gfn’ti,i-u.''n'nuT 
A U )C ( 1 ) =XF0bK ( T ) + I 0BT( T J »SFGS ro j 
rFa8t.T,G4,EG.G)C.n tu 50.) 

:>:i 553 0 = 1,GNTS1Z 
T k'OM5S( i,N) = Tutrr(iv) 

C.iOi'Tu'Je; 

rFC iSt:,Le5,Ey,otGn tu 500 
'tj 5 47 K=l,rtNi 
(FT rPACl,K) = A!;:PiBO(K) 
rQ'-y.i'i>‘>uh 
" I’lTliiffE 

TFtiSF!.,E4.Ey,0)GO TU 560 

.jrtTTEt 1 ,B 33 ) 

P jOMAXCBx, ^THE various ncCUPANCY RATES FOR THE SEGMENT') 

I 'H^5X,^X'** + *»^**»*t‘*t‘’t‘**^***»4**'t‘*’*‘^'*k**t*t»* **»*♦<:♦*/) 

,?4rFUU»2*»6NTSTZ 

•MrtECl ,541) (CQUN'r(n,T = t ,0HTFOu> 

■ FOMuATC 40X,F1 2.6/) 

7 JNTT'-a!IE 
OE'^URrt 

TxD 


Toe routines that follow naineiv o.iTSPi , y*7f S h3 oaa 

are tna optimum ouantiser designs of aaa? and . 1 1 , f ur 

laout o£ uaplacian density. All ot tneoi nave tne sa^e )r?.nt 

Bur4nuT?UE**c‘7 fSRt(Oa*ri«,0(^X0(lT) 

Tnts routine is for ? hit pantisation 
nhUTjOHTXN-lnput to quantiser 
o.iTpUTs U'7‘lT)U’f*0atpnt of tne quantiser 



s 


i'?* 


J 


•’ i 
'' 1 1 

. . j -> 

^ i ^ 


4 i ‘ ' 

4 1 1 
4 1? 

4i '1 

U4 
I I's 
U 6 

41 ? 


r-’f / 1 1 T . * ■' * •» f u. ‘ 1'^:. 

0 iTiHIis-t .510 
?.l M in 
'] 'I'rij'jTsa.IS^B 
‘ ^ 'j I ' m 
? 'T.j'nsi.bio 
'•j n 4ii 

2 •! ';ui11--0,45?fj 

-.. J ‘ ’ 1 T u U ti* 
t R u 
F 


T'j 4Ai) 

)Gn Til 411 

AnD.f u»'rTU.GT.f-.9Bl6jm'l fn m 


S . ) k n b r I N t b- N i; s H 3 (. 0 M' I , 0 iv.T □ 0 1 ) 

is auantlsatton 

Tfj cbT 

] f ^ J \ R ', ' J i ' T , 0 , U 2 . A N 0 . ( y M T T i '.| . G p . C - . b 0 4 ) ) ) G a T □ 414 

1 1 r R M I ■ '2 • “ • ^ ^ ^ ''* « ^ • ’ • ’ i *> > i ^ 

,J, I pj ^ 1 !! . fj' T . ( - 1 . 1 H 1) ) . A r ’ D . U) 'i T 1 w . G fe: . t - 2 . n 5 m G n T 0 4 1 b 

Gj 'i'n 4 17* 

') nLi!!Tsii,222 
Gj jn 111 
0-iTaUTa(i.785 
GJ TO 41? 

O.raKJT^l.bTb 
GJ Til 417 
O'vi "o'V *■=? .994 
GJ XO 417 

9iM’0nT=-,222 
Gj TO 417 
0'jTQUT3:».78b 
GJ TO 417 
0',fDilTs-l,57fi 
CJJTTNGI;; 


R^iD 


;Ti/Ri^ 


SU«RnilTIMB yl«TSR4{Of<TIN,ORTtJffT3 
Tjls roatiae is for 4blt quantisation 
T F f ( 0 N f 1 H , G E . i> . 0 ) . A r • r . f U H X‘ I H . L T . 0 . 2 6 6 1 ) G 0 T U 4 .1 U 
If E ( ( OiM t i »< , Ge , 0 , 2^6 ) , A GO . C yMflil , GT . 0 , 566 ) 1 GO TD 4 1 1 
ir( (OhTI?1.Ge.O,56t»}.R;..D,Ca?4TTM,u'P.O.9iO))G0 TO 4l2 
TEC t0RTlW.G£,a.910),Al<D,(qnTa.MT.1 ,33:7))GO TO 41.3 
1E( lOJTI^'i.GE. I .il7).Ah<D,fatfi‘l;i.l,T,1.821 j )GQ TO 4l4 
r E ( ( OMTI , GE , 1 . 82 U . A wf) , ( yMt IiJ . uT , 2 , 499 ) 1 GO TO 4 1 5 
t K C ( ONTl . G£ . 2 , 499 ) . AwD . C Q« ft M , tiT , 3 . 60 5 ) } GD TO 4 1 fo 


I 



1 ? 1 


■1 A 1 

1 1 ^ 
1 i 4 

1 iP. 

'M7 

1 ,) ' < 

1 I 

1/:^^ 

1^1 

'i2 7 

■r'/7 

rii 

12S 


r f f '.} '! l T U . G R , 3 , ft u B J r.o I’f) « 1 7 

f r f A t i * l-'T • * 'U * 'i"'" :• -GR . 

r ' f !' I . * t ^ S’ ) '' . ^ '-i n , ( g T T f. 

T-f J f ‘ •I''*' t-.Bfoft) 1 . AND, fyNTTr.,GF:. ( 
Tf f (.0>.1 1 'i.ril*. l-.'liO) ) .AijD, fv>MTT -J.GF. ( 


(■-,26(0 )1.,3 T,, 41 « 
.GF. JlGD Tu 


4 1 ^ 

,91 :J1))G0 Tu i?U 
1.3 17JUGJ TH 421 


T »- f f -Si -T’ r •(? r ,« > r 4 ,4 i V * 1 I. 4 , vjr. , A - i . 5 i f j i j J 1 m 


IRC (.71 Tl ) 

( Oi-tM 
7 .rT-J JT=«»4 

J I t s 

T3 4 25 

2:''"fUt'4=:0,40? 
Gj n 4 2B 
7 .rjnT=U,7?ft 
'^J I'T 4 2B 
7 aTJIITsI ,U9;i 
7J J'7 4 25 
2- '’J!M'si,.54o 
7J rn 425 

70 TO 425 
0 A 1'J'I1=:2 .8 05 
o.j X^'I 425 

O'J Tjii'i S4 , 3 1 o 

7 i 10 425 
0.,TjifTt=<-.U& 

1 10 425 
n'Jiris-,407 
73 TO 425 
ONTonis. 72e 

73 TO 425 
OTTOUTs-i.095 
70 TO 425 
O'iTUilTs*! .540 
73 TO 425 
Q^TjU1=-2.103 
7u TO 425 
03TQt1Ts-2,895 
CU^TTMIK 
RSTURft 

FND 


<i . 1 - 1 . 8 7 n J . A»U! . I 0.; T . Gfi 
r.i. (-2,409) j .A^'u. c owTio.rtfc;, 
3 lb 


(-2. lOy) ) )(;il TU 
(-3,a05) 1 )G0 r,, 


4 7 2 
4 7.1 
4 7 1 


0 


BUBRnUTlME l)ECODE(nECXN,M,X,MAXQPD,DECDUT,SEGST0,IKfc;t.tl) 
Tnls routine DECODE accepts tne deouantlseiil sa»Die value 
and reconstructs the final output based on the or.evious 

0 U t P il t S 

IMPuTSs DECI«^:-lDPUt to tne' .decoder , from tne. dequantiser 
■'1-Tinie no.of sarules In the seqment 

x-Vector of tan oalns of the linear oreaictlon ffirer 
"^AXOPD-l -Faxliiiuiri uRD-ler of tne aoove filter 
5«,”aTD-5tandard deviation of the segment 
iBRuRi-Seiection Keyed in by the operator for adapting 

0 JfPiJl’j DECOUT-Final analog output of the decoder 
OlO£^St07 0ECi«C!00) ,,K(2O>,UECi0aTCt00) 
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[>' f I 1 , 0 ) sFUSTHsI .'I 

}j .i. = l,.»! 

( f. .^0, 1 )GU TO S35 
', J TO *>40 
0.4 0 ''.0 : CT.'^Oc: 

0 ■ : C U " ’( ( 1 ) = ' > Si C £ 0 C I ) * 6 F U s r 0 

G j .1,0 5 40 
OO’’ 7 00.CTw:(£ 

i 1 f'.>BK = 0,0 

oj 55a J = 1 , ti'IAXdRU-l > 

!;FC([-0) ti6.0jr,u TO 555 

'U'^OjsXC J)*OECunt(l-J> 

Opr)i3K;;!)KDBK + nXFl> 

550 GO^lTi.-Iii: 

555 C0OTi:iOIi 

OiiC JO T CDs ( UFC loti) ♦ScJCSO’U + JFOHn 1 
5,jO CjMTIi'-Jllb; 

» FC.l 5hl,(i:i .Eo,U)Gn TU 560 
T O J 560 ,£ = 1.0 

? oexaoTC I )=6EConT( tj^st-xs'^o 

5tjO fOOTf (.Mfh; 

F;in 


6i 0 
615 
620 
625 

630 


SUBKOUTlMe CnUHX4(0oTUr’I,CnuMT) 

rnls routine C 0 UfiT 4 oives tne quantiser occupancy statistics 
for the 4 olt quantiser 
I jPUT:OWTDUi’-OutPut ot the quantiser 

nuTPUX: cnuNT-Vector contaminq the quantiser occupancy 
statistics for the entire interval 
nioxoTO' COUNTClfe) 

; 0,126)00 TD 610 

TFf‘.o iO.iT.Ey, 0.407) GO TO 615 
IFrU:'i’10UT.EO,,r26)GO TO 620 
TFCQJTDifT.EO. 1,095) GO TQ 625 
TrfONTOUT,FQ.l.5^0)GO TU 630 
rF(0M'ihUT,Ea.2.i03)GO TO 635 
IF(yMTOUT,EQ,2.895)Gn TO 6^0 
I F ( y 51 ’ 'UT. F0 , 1 , 316 )GO TO t:>Ju 

•: i J'.ii .Ea,C-. 126 ))Gn '"u 650 
'>G' 0 '%EQ:t-. 407 ))GO 6? 5 


' iCij'r ,EQ.£»,726))Gn Tu 060 
' i.'ir',F0,(-UO95)}GO in 6b5 

fT.tr,KO,t-2flOi)5GO in 675 


TF( 

IFCv 
TF('.; 

IFf V 
f... 

TFC^irrouTiESIi-^iegsiiGo to 6bo 

TFCQMTOUT.EO. 1*4.316) )G0 TO 685 
CjU«TU J5=Cnoy’T{U+1.0 
Gl) tn 690 

Canirn?)=cnuTT<' 2 ) + i .0 

r.) TO 690 , „ 

C yf frt T 1 3 ) SCOU h t C 3 ) + 1 , 0 
GO TO 690 

GjOhiT( 4 )=CnuhTt 4 ) + l ,0 
GD TO 690 

CO*GFr(5)»COliMT£S)tl ,0 



<4 

G J I 7 6 9 n 

C )■' /rt 5 )sCn(;‘!T(t>) + 1 .h 

'•J i:i 590 • 

IS i 

7 J r , !’( 7 i=CQU 0 I'C?) + 1 .0 
• J rl 590 

' • J I t 590 


” j".>/r(, 9 i=conN'T(y) + i .0 

GJ in 690 

7 J'):iT(lu) = 7 n!iwT( } o) + 1 .0 
73 TO 690 


fy t 0 

7J0..I f'cu laCUduTCl i) + i ,n 
70 l-< 590 

fMjS 

7 jO;.T(i 2 )= 7 unf..T( 1 i') + l ,'1 
7 } iO 590 

IS / (> 

7 JO.rrci j)=ntiU.vr(i 3 ) 4 i .0 

G ,j '! 0 n 9 0 

f ; / 

7 .jU'T'’( 1 4 laCO'Ii'iTtl 4 )t 1 .0 
73 TO 590 

6 H '' 

CJi’rtTUb)=CUd0T(15) + 1.0 
7 J TO 590 

^ r% 

ron,^ I ‘(1 6 )=Cijn .ni 6 ) + i ,0 

fS iyr '' I 

7..0Tru'IE 

'«h/l’u!7v 


‘JUfliinuTiMfe h>JCnD4cni 
t'nls routine tiWcniH 
HufitT'ian encodinq 
r 'iPuT;0i'jT0UT-0utPut 
GUTpUi'SjiQUT-ls tne 
output 

TtNTRU-is tne variable lencrtP coded vector of tne 

t)l'n3iTy' codi€J 

MNI-Slze of tPe variable lengtn code vector 


.;TUnT, inuT , Tfc;f^TflU » «Hi ) 

does the bluarv and variable ienoth 

0 "fc t* 0 Q u c 3 i ^ tt i s 0 r 

binary encoded vector of the <judPtlser 


Pi men ST OH 
TE(0^JiOUT 
re" (O'." TOUT 
If (QT TOUT 
Tf (OHTOUT 
Tf CyUTOUT, . 
Tf (d'n’OUT.EO 
Tf 

tf Cvv i''^U”'.Ea 

ioir^-,E0 

TFfQ^IOUT.EU 
!f CO-’TOUT.E'^ 
Tf fO -rmUT.EO 
rf cy.i roiiT.ey 


rOUTf4),TE*' 

E^,0,i?6‘)qC 


Hipocn 


.U,40/}GO 

,0.??6)GO 


E 

|o;i;o9siGn 
EQ.l .540)GD 
” 2.103)GO 

2.895)G0 
4*H6)Gn 
,1263)GO 
:4075)GO 
,7263)Gn 
1,0961 )GU 
1-1 .54i))3GU 


'•i 

TO 

TO 

TO 

TO 

TO 

TO 

TO 


710 

715 

720 

725 

730 

735 

740 

745 

TO 750 
TO 755 
i) 750 
10 765 


rFfy^lTdiJT.EO, (-2,103)300 
" " .(-2,e95))GU 


TO 

TO 

TO 

TO 


tf fOHiOiiT.EOv , , . ... 

TF ((?^*TOUT. Eg, (-4,316) )G0 
71 r tjJM'dlal 

t:o'ti2)=o 

TJHTOJsO 

Id0'i'f4)»0 

C kt . G V A 0 ^4 4 0 1 1 00 1 f 1 Er^TKO , Km T 3 


770 

775 

7bO 

785 



17 '* 


7 1 




7 i 


7 3'-. 


7'r 


7 4S 


7 51'* 


755 


760 


lEKTftn/'wT j 


CJ li) 7y{) 

I J!'l f i) = i 
' in ( f;^)=t) 

( JHi f n=:o 

V , j i I j (■ ft ) s 1 

C Af.l, vAL'i4 1 C tout, iRriTKn.^n'.T j 
'i I i ') 7 y 0 
T J''Tfn = l 

T )'!Tf2)=0 

T Jii-i f y)=n 
C.'.r.L. V l TGl!T, 

'J Tn 7 90 ' ' 

f jnK Dsi 

r:jni(2) = u 
i. j'ff (i) = l 

r Jt’iC4) = i 

^ '1] f ’ f|-, ^ L J ' 4 3 C 1 0 U T , I F i'i T H f ' , r , * T 1 

'l’jnT(U = i 

T 071 ( 2 )al 

I JHTf i)=0 
T JiiT(4) =0 

r a{,\. V a 6M4 4 ( j uUT , IFhTHTl , r‘i. T j 

(7.1 I'll 7 9(1 
TjilTd ) = 1 
r Ji)Tf 2 ) = l 
r j*i'i (31=0 

r :)nx(4?*i 

CAf.t VAG>!45( TOIfXMFNTKO,! :.)T) 

Gj to 790 

raoTf nai 

7001(2 5 =J 

T3nT(33=J 

TJ!!T(4 3=0 

CAC.b VALN46(TOnT,IEivTRO,M;vT J 

Go TO 790 

IJOK i)=i 

TJOI (21=1 

TJUT(3)*i 

TaoTr 4 )*i, 

CAI.L VA6W47(TOUT,iEwTRrt,Miv(T) 
G:j to 790 
TQ!(TC1)=0 
raHT( 2 )=i 
I:;)itxC3)sl 
raoT(4)si 
CAr.f. v|gN48(IOUT 

r J ' < ( n«6 

roox(2)=i 

'‘■0OT(3)sl, 

r,)»,iT(4 )=o 

C '\?.b V A lj?i49 ( TOUT , 1,F> TrO , '5vT J 
GO TO 790 
raorc i )»f) 

7 .J!'(TC 2)«1 
rofJT (3 3 so 

( jiiT(4)sl 

C i\ 5 V A li I ‘i 4 A 1 1 U H T , I E T R 0 , M h i j 

G,9 to 790 


lEhTHn,MbT ) 



f it :t iC - 5 r 




T jn i ( ^)=1 
i j"! {3)=n 
I 3iIT(^n=0 

C Af.b V TijUX, ) 

'l.i I'O 7 30 
7/' fji'.cCjiisn 

f -M' i f/i)=0 
I.iMI, C J)=1 
r );IJ (4)=1 

7 a r , L V "a b f ' 4 C C T u 1 1 r , 1 F , . T ; < n , ?■' . -i T ) 
aj i':l 7 30 

7 r". r )=o 

r 2 ) =0 

i: jni(3) = i 
TJOf (^)sO 

C.hIjL V ALf'i 4 0 1 1 u!f i' , 1 F.V THH , F , T J 
r;o TD 7 30 
7h'‘i (oU'i'(l)=0 

rj*'i (2)5=0 
T jd'i (3)550 

T J'lT(4)55t 

r n f , I , V A IF’’ 4 K ( T 1 1 (! T , 1 F i . T K n , ! ■ A' I j 
0 J TO 730 
Ir.'^ T,|'tr(l)=0 

T'jOi f^)zO 

T,)nT(3)=0 
I jllT( -4)s=0 

C Atjb VAUW4F(TDUT,IFN'rHn,MoT ) 
7'F.' COMTIOUE 

R E T U R W 

Fi'jf) 




R i5 


8 20 


SUBHOUfUJE UECFl4(ITRANs,M,aECI.|) ^ 

Tnis routine DECFI4 is tOe decoder whlcn. tiaps tne binary 
code, into tne respective dequantised samples 
T'iPUTS: ITRANS-Binarv encoded vector expressed In 2 

dimension for tne segment under consl- 
f 1 0 :F iS ti JL 0 n 

M-Mo.of samples in tne segment . 

OUTPOTt DECT^i-’Blnarv decoded sample that is input to 
the decoder . , 

DIMLMSIUN IXRAfJS( i00,4) ,DeCrM(luO) 

no BOO 15S1»M ^ „ 

If f ttranr(t,u,eo.o)c;o to fis 

TFCi ■aA,.Su4i.i0.0|r,O TO 820 


*U r,3 

F(r,4 


316 


T F ( I 5 'r' « 

TFdi'HA 

neciMCt js4 
c ,0 to 890 
TFf lfRANS(T,2) 
!FCiTHAMSa,35. 
TFf ITHAi’^S(I,4) , 
f:)£CrKI)s-0.126 

GO to 890 

TFf ITRArlSCl ,3 

TFCri'aAi'<S(t,4 


S6;0)GO TO 825 
£0.0)00 to 830 


EO.OGO TO 835 

EO, 0)00 TO 840 
eO.O)GO TO 845 

EP. n)00 to 8S0 
EO.OJGU to 8o5 




l?o 


-J j ■) 

q/jS 

I) ’) » I 

I? 1,1 1 


U 


Q / ' 

;s 

HbO 
9 9'’ 


ifcO.ojr.Li in 


.to.njnu 

. t.0 . 0 ) j 


TO 

;io 


OiiClMT ) = 1.095 
TJ rn 990 
tfd rrtAAi,9CT,4) 

1) n 990 

r)=2.H95 

•7,J iO 9 90 

If (lTi<AiviS(T,3) 

If ( I <'KA,jS{ T 4) 

ofo I n r)=-t .540 

‘d IT 890 

Tf riTHA‘l5(T, 4) ,H;o..n)r;ij 
')';ctu;rj=-o'72§ 

Oj 10 890 
0£Ci:'ItI)=-0.407 

c;j n 8 90 

Tf (l rKA4S(T,4J .b;0.O)OO 
oyrpnt j=a.40/ 
r,J TO 890 

ogriotDso.??!! 

04 TJ 890 
) t’j r 1,9 1 1 J s: 1 ,540 
9J TO 890 

T f ( i f H A 4 R n , 4 } . fc;o .0)0; d i'p 
Ofri 4(T)=-7.995 
0,) T'1 8 90 
OliOI -itDs-^.Kd 

0) TO ft. ^,4 

oari .u us-i .095 
0,3 xn 8 90 
niiCiMc 17 = 0 . 12 b 
nu xn 890 
neri:im)*-4,3ifi 
C 30X1 HUE 
CO^’TTMME 
RETURN 
Ron 


«rt5 
8 ?0 


rn 87 5 


in 8 8 0 


8H5 


Ttie follot^Klno lb routines are tne variaole ienoth 
coder routines nam^d as vaLn 40 to VALiHr^tOe last 
letter oeln'.j i.n 'oevadecifnal ot u to Ib.fhev have the 
sane Input output ocirar'eters ,T'Oe "4"in toe. naRteof tne 
routines indicate tnat tbev belong to 4olt yuantiser 
SUBROUTiUE \/AEW40aonT,IFNTH0,MNT) _ , ^ .... 

TmPiiTj TOU un'v encoded vector of tne quantiser output 
'-i; li'Ui-v'ai'lable lengtn coded vector of toe olnarv 
c o d. ^ 

MSi-size of the variable lendth code vector 

OL^EUSTOfl I0U1'C41 ,TE«XROe75 

I'.f'-'-.iu-'uiDcl 

i.f 

■'^,,1:^7 

, , , , 




! 

1 

1 

i 

I 

127 i 

j 

! 

.'^iHq nirn ke; vM/f’-ti t ujt't, j 

'r ^ I’ c 1 ' ’ ■ T M > , T fa: '1 K U ( 7 ) 

r t!.w 1 W;j ( i )sf) 

Te;:''TWfj(k) = 1 

TB'ilROf 

1 '1 { 4 1 =1 

'Vm' 1 S.1 
“■ T 1 ; u , j 

p:,p 

1 

! 

\ 

SUBKHLlTlNt; VALr''42t TUB'j>IFrvTHn,‘’,M ) 

TOUT( .-n ,Tfc^’i Ru( 7> 

T 11, ' 1 P u ( 1 ) = 1 

T£MTRy(i')=0 

ieMTPOCi)*! 

Ti:;«'lRu(<l)=0 
*'1 is T s 4 

Rt. Tun:^ 


S U B R 0 u T 1 « t V ft L f ' 4 3 c T Q ! J f , X F N T K D , M T ) 
OX'flfalNSiaN IOU1C4) ,IEF'XR(J(7) 

I£f'!TRu(l) = l 

TSMTROC2)=0 

ts;’7XR0(3)sfl 

TS'''ilROC4) = l 

Itt;^lTR0(5)=0 

''l.'MisS 

MfiTURrM 

E‘'JD 


S 1 m H 0 U T X « E V A Ij H 4 4 ( I u U T , I F N T R a , M iO J 
DIMEMSTOS I0U1(4) ,IE^3TP□^) 

T£-m,pG(1)s 1 
ri<vViRijC2)so 
ieNX'RO(3) = 0 

T£?^TR0C4) = l 

T. 3(5) = 1 

'i '■ f s S 

HEfURW 

E40 


n J 3 » 0 U T I E yf A L 54 4 5 ( T 0 0 T , I x: rt T H a , «s ) 

D £ M E ,'l S 1 0 M T C i ' J f ( 4 ) , 1 fa. • 1 P 0 f 7 ) 
t£rrRDCi)=i 

tSfFfR3(2)»r> 

T-2WTiRD(3)»1. 

1 

f 

1 

1 

1 

[ 

i 

1 



]■ »lU 5)=.1 
T.i'' I’!' j(b'»=r) 
1 lT = r^ 

i’’ .1 ! > 


.1 '! M i ' ' y ^V'- V * I-' ^ c T u u r , I F .M T .o , M ,v n 
TUUi(4l 

I f:, 'i r '< U ( I ) = 0 

rc:-'i.v:?o(2) = i 
rirf 7 Q( 3 )=o 
Ts u'^ac i)=o 
T£’>!l'^>U(5) = J 
T£vr^j(fj)=a 
' 4 U = 5 
7fiTuRM 
^ V) 


SuBKlUTlMc; VAL\»47 (Tunr, lEw’fm,MwT ) 

'HMtisiiJj rooT( 4 ) ,TeMrRaC 7 ) 

T.e''TRaci)=i 

T£''iXR0(2)=0 

t£4TR0(3)=t 

Tt!:‘'T«0C4)=:l 

iavTR 0 ( 51 al 

lE'viTROCS)®! 

i** t s 6 

RETURN 

Ei'ID 


SUBROUTI^lE ;/ALM 4 RCT;::' 

Ol^ifc;MSIOA! I0{JTC4) ,lii 

T£KXR0CU*0 

l£.MTRO( 2 )sO 

MX 1*2 

RETURN 

EiXD 


t, r^-.:THO,VMl) 

■ipor?) 


S u B H 1 Li T 1 7 E ¥ R L M 4 9 ( 1 0 n T , I F N- T R 1 , M X T ) 

DlMfe;xsT07 It''tJtC4) , TfjPVRiHR) 

T£f’''i‘RO(l)«Q 

TENXR 0 ( 2)*1 

fgNXRDCi)st. 

fV.ilsR 

returu 



SU^h'U'Tl-vt-; V lOI'l j 

” i. ’^K ! 3 I ItiUT( 4 ') , Rli( 7 ) 
rj'/i'rs'ri ( i) = i 
1 RUC 2 ) = 0 
IP rj( j) = o 
T^r'TWuf 4 )sC» 

■■^jTr 4 

!•■; 'j n 


5 u !> P n li T i N f V A L M 4 B 1 1 om , i F <4 T K n , M w T ) 

rnf'’li'!STOr! IUI!T(4) , Iff'TBOf 7) 

TifilRUd )=s 0 

T6;«TRU(2)=:1 

Ife:''iTRU(i)=0 

f t;'iTRC)(4)=f> 

Tc;’'TBU( 5)=0 

RST’UP.i. 

FivO 


S l m k n 1 1 T I M E V A L M 4 C { 1 0 U T , I F N T R 0 , M k X ) 
ni'iEJ'JSIUN IOUTC4) , IEWXPU(7) 
TENXROCnsl 

IEMXR0C2)*0 

IEMTPO(3)=1 

T:V- 7 RtiU)«l 

TE,"} 

- ■('RL«( 6 )aO 

' ,, t -S 


SUBhOUTI^E V AEfj 4 Dt| 0 in#lFf|rR 0 rM«T) 

D I me ‘33 r OM lou r C 4 ) , TL* ( 1 ) 

TEMlROd)*! 

I5NTRO(2)»0 

Tei^TRa(|)al 

lENIRoUl*! 

tE'MPuCSlsO' 

TE l.'kljf ! 3)=1 

R. 



I 




1 ■,! i 


O'-jO 






'US 


QoS 


9/0 


0/S 


980 


' ' 0 i la( = t 

f U. :;Ui:3)*:i 
I ■’*i.l;U 1 ., •■(:) = 1 

0 i I'D Qy.) 

UU lF.lTntAtT ,4 ) 
t !i5U,l.) = l 

. On JSf i,21=0 

1 [’'U\0bC 1 , i) = 0 
*■ •' "j 3 ( 1 , 4 ■) = 1 
' - 1 n q 'j 0 

u’l iFj, ri<M.r,5j 
^ F c t F i r 1 < A ( T , s ) 
UUiS'xTrtAd;?) 
T I’^A-lSC 1 ,l) = o 
r I'OADSCl/iDs') 
f roAssci,,i)so 
n'0A>(3( It'D®! 
li TH 990 
I U>A''U,(1,1)=0 
>OA:'lftCl,2) = l 
f 'nu\'l3f J ,3) = 0 
UqA''J«( I.4)®0 

f!} I'D 990 
T rSAfUU L, j ) = l 
f 1‘MAM.SC I ,2) = t 
’ rPAHSd, J) = 1 
T i:OADS(l,4)B0 
9J in 990 
Tr9A93(i,n=a 
r 1’9 AiMS( I ,2 )=0 
TrPAM6Ci,3)=a 
rxRAMsr If 4)30 
so I'D 990 
TFf IFlTRAd, 3) 
TF(.lFXTRACr,4) 
TFdFlTHA(T,S) 
TF(IF1TRA(T,6) 
H'RAMaUf i)3l 
T TRAMSC 1 ,2)st 
TrRA43a,3)sl 
Tr RAMS (1,4)31 
'.0 TO 990 
TFfiFITRA(;i,4) 

i:fcxf:itha(t, 5) 
TUUMSCI, U»t 
T TRAMS? 1,2)31 
T trams a, 3) so 

rTRAMS?I,4)sO 
GO TO 990 
TTRAMSCI,i)*t 

r r»A:Ms(x,2)«o 

ITRAMS(I,3)sJ 
r lRA^'SCl, 4)»0 
GO TO 990 
ITRAMsaajsO 
I '3 ( 1,2) «! 

! r=-£-.DsfI,3)30 
lfRA«S(l,4)s1,^ 
GO Til 9 9'- 
I li‘RA?U5(l,'Usl 
ITHaMSC 1,2)30 


GO. ft) GO rn 04 0 


. t '^ 0 . 0 / G I ) 

.tO.O jGu 

.t,0,n )C.J 


fO Qh:> 
IP 9SS 
TO R50 


. GO . 0 ) Gu 
.GO.O)GO 
,eo.o)Go 
. GO . 0 ) GO 


DiO.OJGU 

go.oJgo 


TO 0o5 
TO 97 0 
TO 985 

rn 999 


TO 975 
fO 98 0 


m 




0 ■) 


1 < > 1 u 
1 (- 1 
1 0 ? 
KJ?b 

1 0 y 

1 3 :> 

1 1»4 a 

\ 04i 
1090 


T i”Aiy( l,2)=0 
T fPA'isf r , 3) = 1 
r rHA^!b(i,.i) = .i 
o ,( 'ni oy6 
T TKAMiSd , l)~1 
'f 

i I, j)=9 

r r'-f/wsc 1 , 4 ) = ! 

rj^'f i ivtjK 
" J '■! T [ 4 U fc; 

!•: > 


oc cup Stacy 


s ..t n H 0 u ’’’ 1 N fc, c n in ' t 3 1 0 !.■ t u i it , c n i j m i ) 

'I tils routine cnu^T3 nets the Muantlser 
^tj’tlstxcs lor the ibit ouaattser 

the quantiser 

'u • I'Ui • ^„fi jfix- V ect or containing the quantiser 
^ statistics 

blf’bMblON COiit'jT.iii) 


occiipancy 




■n^y- 


o.:i'72)(.in 
0,7P&3G0 
1 .‘j7t,)Gn 
2 . 994) GO 
(-,222pGO 

c-.nsnGo 


t-1 

(- 2 , 


576) J GO 
994)) Go 


0 


Tf fQMXOUT.Fv 
fFfyMlDUT.PO 
TFlQM'iOUT.e 
TFCUMXOIIT 
TFigMiOUT 
TF(0 ITOUT 
TFCyPTOUT 
TFfVMTOUT. 

Gon;uHi)=cnu9 T(i) + i ,0 
GO TO t09U 
C0Ui4Tl2J»CDUtJT(2) + 1 
GO TO JU90 
Crj!l^iTt3 JaCOUNTC3) + 1 
GO TO i 69'0 
C0nqri43sCncjNTC4) + l 
GO TO 1090 

ronfjTl5)rCriUMTr5)+J .0 

GO TO 1090 

C0l>MTt6)8:C0Uf'!T(6)f 1 ,0 
GO TO 1090 

COI'WT(7)sCnuMT(7) + J .0 
GO TO 1090 

Cj(liMTl8)aCOUWT(8) + l .0 
roMTiNue 

0!;;TUR^ 


TO 

Tl) 

TO 

TO 


1 0 1 il 
1015 
1020 
1025 
TO 1030 
TO I0i5 


in 

TH 


1040 
1 045 


,0 

.0 



11 1 V' 

U 1 i 

I 1 7 

I I i 

1 llu 

1 1 

1 1 4b 

1 1'4U 


il!§ 

jnur,&M.2,,994)Gn to 112 *; 

r p f i,r ! 1 U J 1 , fi V , I » , 7 b ) (.n ili*^ 

T H* f r ■^, .'4 4‘ n I . T 1/ ^ ^ T y i V ^ - I i 


in 

■| n 


[ - J H n ^ “ I • ^ v> : ( - 1 . b 7 0 ) y G u 
(-2.yQ4) )ru 

r .j'ii (2) =o 

r Ji' i ( i)z;) 

" ' ■ !/ ,, ^ 9 3 0 { T 0(1 T , 1 F>» T kO , h ,V r } 

'1 i ! 119 0 
TJ"1 ( 1 )S1 
T Jt'i(?)=0 
fjf'i(j) = 1 

1 ibiii,iFi,THn,f’NT 1 

GJ J'O ll9u 
r<jnT( 1 ) = 1 
1 011112 )al 
1 JUl(i)=0 

CAT.L VAL.:!,l?t rU('T,lF(:THn,M.',Tl 
G.J fO U9u ' 

T jnxd )=i 
nnif2) = l 
T Jill (3)=1 

c AGL^V JL^Ji ? ( tout , IFmTKO , Mw T ) 

tjnx(2)s} 

ionT(3)=i 

CAf.U VALli34t TUHT, lFlvTHn,f»!,T) 

GJ TO tl9U 

I JIlxCl )=0 

T0t!T^2)al 

IJI1T(3)=0 

CALL V A L r< 3 5 ( T 0 1 1 T , 1 F pj T H n , F ,4 T ) 
GO TO U9U 

roiiTC 1 )=o 

I JHTf 2)«0 
IJllTCi)*! 

CALL. VALM36 ( TOlUMFilTRn, Mill ) 

GO TO ti90 

lOUTCl )aO 

r OUT (2) SO 

TJUT(3)*0 

cor.b \/AuM37( lOiJTf IFN'l’MO,’A.T) 
COMTtNllii: 

RETURN 

SNO 


I Hu 

I I 4b 


fue foUot»?laq 0 ro'ittn'^s are t!Ti«^vjir laoie lenatn eo<S€r 
routines natnea as vat.mSu to /AfdlJt.fne “S'* in their 





refer to jhit nuaot 1 ser . 

* ‘"output vector of toe quantiser 

^ n p Ji >: [k,'I I'RD- vai'i^bie lenytti encoden vector o£ the 
oioarv Code 

S./r:Kn,Tl.’i;' vector 

r i'Ou( ^)=0 


'< T s ? 

dt,T(jpr^ 

K . o 


;*> U H !• q u T 1 fi t: V ^ LN 3 I I T uirx - i F l^! TK q , M T 1 

JUnTfi1,TttK'0f /) 

I b ' 1 R u V I J s } 


Tb''iOD(2): 
rb ’iRu(4): 

U;vTs:4 ■ 


■S U p. F q U T i N b V & L 3 ? (, T LJ tt r , t P .'J T H D , « .-vi T ) 
'UMtdSTQ'V f U'JTf 3) » Tfc;^'TRUf 71 
Tb"brRU( i ) = 1 
tbM iRU(2) = t 
rSM TRij(3)a! 

ToUi'R0('i) = 1 
T!:;WTPa(S)=0 

RgfURM 
K ;< D 


S U (3 K D U T 1 t'f E V ft LtJ 3 3 t T UUT , 1 F N T RO , M i'l T ) 
DIPEWSIOM IUUT(3) fTENTRorVI 
l)a1 

^ f.i ».’ p f 2 ) ® i 

'T!m.,K4;»i 
r,:"VPij(53«l 
1 :';:''i;‘*'i;Kb)xO 


l«6 

RETURN 

EMO 


.SaOHOUTi^^E ^ ,niaM44 ( TOUT # IFoTRO , Md T ) 



C M'M M M 1 


1 ^ i 


'.) L ’p, I u . 
I ^ > f'K'JC I ) 
M , T = 1 

V . r , 


i) 


v/ AljfLiS ( TU^' J', IF'/Vivn 
|0nT(3) , TENlPOf 7) 

1 b. ' I KjC 1 ) = 1 
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